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ABSTRACT 

We present the results of a comprehensive Spitzer survey of 69 radio galax- 
> ; ies across 1 < z < 5.2. Using IRAC (3.6 - 8.0 /zm), IRS (16 /zm) and MIPS 

^ I (24 — 160 /im) imaging, we decompose the rest-frame optical to infrared spec- 

Csl ■ tral energy distributions into stellar, AGN, and dust components and determine 

Q ' the contribution of host galaxy stellar emission at rest-frame H— band. Stel- 



lar masses derived from rest-frame near-IR data, where AGN and young star 
contributions are minimized, are significantly more reliable than those derived 
Ph' from rest-frame optical and UV data. We find that the fraction of emitted light 

O ' at rest-frame H— band from stars is > 60% for ~ 75% the high redshift radio 
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galaxies. As expected from unified models of AGN, the stellar fraction of the rest- 
frame if— band luminosity has no correlation with redshift, radio luminosity, or 
rest-frame mid-IR (5 /xm) luminosity. Additionally, while the stellar H— band 
luminosity does not vary with stellar fraction, the total if— band luminosity 
anti-correlates with the stellar fraction as would be expected if the underlying 
hosts of these radio galaxies comprise a homogeneous population. The resultant 
stellar luminosities imply stellar masses of io^^~^^-^ Mq even at the highest red- 
shifts. Powerful radio galaxies tend to lie in a similar region of mid-IR color-color 
space as unobscured AGN, despite the stellar contribution to their mid-IR SEDs 
at shorter-wavelengths. The mid-IR luminosities alone classify most HzRGs as 
LIRGs or ULIRGs with even higher total-IR luminosities. As expected, these 
exceptionally high mid-IR luminosities are consistent with an obscured, highly- 
accreting AGN. We find a weak correlation of stellar mass with radio luminosity. 

Subject headings: galaxies: active — galaxies: high-redshift — galaxies: evolu- 
tion 



Introduction 



Luminous radio galaxies were the first class of obscured, or type 2, quasars to be dis- 
covered and characterized. They have accreting super-massive black holes whose continuum 
emission at UV, optical, and soft X-ray energies is absorbed by dust, thus allowing a clear 
view of the host galaxy. More recently, hard X-ray and mid-IR surveys have identified the 



radio-quiet cousins to luminous radio galaxies (e . g.. iNorman et al.l l2002l : IStern et al.l 12002 



Martinez-Sansigre et al.ll2005l : IPoUetta et al.ll2006l : iLacy et al.ll2006l ) . The main evidence that 
radio galaxies host super-massive black holes comes from the high luminosities of their radio 
lobes, which are fed by rad io jets originat i ng at the host galactic nuclei. The lobe spatial 
extents (up to a few Mpc, ISaripalli et al.l l2005l ) and luminosities (Li.4ghz > lO^^WHz"^) 



clearly rule out a stellar origin for their e nergetics (e.g. 



tation unification scheme for AGN [e.g., iBarthel 111989 



Antonucci 



Reeslll978 ) . In terms of the orien- 



19931 : lUrry fc Padovani 



19951 ) , radio galaxies are radio- loud quasars seen from an angle where an optically-thick torus 
obscures emission from the region closest to the central engine. 

Due to their large radio luminosities, radio galaxies were the predominant way to probe 
the distant universe until the advent of 8 — 10 m class telescopes and the Lyman-break tech- 
nique in the last decade. In fact, ra dio galaxies were the f irst galaxies to be found above 
redshifts one, two, three and four (see lStern fc Spinradlll999l . and references therein). Across 
cosmic time, the host galaxies of powerful radio sources appear to be uniquely robust in- 
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dicators of the most massive galaxies in the universe. At low redshift this result has been 
known since the first optical identifications of extra-galactic radi o sources showed them to 
be associated with massive, giant elliptical (gE and cD) galaxies (IMatthews et al.lll964l ). In 
the more distant universe, indirect evidence that this association remains intact comes from 
the detection of host galaxies with r^/^ law light profiles in Hubble Space Telescope (HST) 



observations of high-redshift radio galaxi e s (HzRGs) at 1 ^ z < 2 (jPentericci et al.l 12001 



Waddington et al. 



2002: Zirm et al. 



resid e in moderately rich (proto-) cluster environments (ILe Fevre et al. 



20031 : iBunker et al.l 120071): the tendency for HzRGs to 



1996; Venemans et al. 2002, 2003, 2004, 2005, 2006 



Stern et al. 



1996; Pascarelle et al. 



200 



^; the spectacular 



(> 100 kpc) luminous Lya haloes seen around several sources, implying large gas reservoirs 
( iReuland et al.l l2003l : IVillar-Martm et al.l |2003| ) ; sub-n am detections of HzRGs, implying vi- 



olent star formation activity up to ~ 1000 M0yr ^ (lArchibald et al. 



200 



2001 



Reuland et al. 



i^ ; and a few direct kiii ematic measurements of HzRGs ( Dey et al. 19961 : Dey &: Spinrad 



19961 : iNesvadba et al.ll2006l ). The most compelling evidence of this association of HzRGs with 



ble, or K — z, dia 


2;ram for powerful radio s 


ources 


fFie. m Lilly & ] 


joneaii 1984: Best et al. 


1998: Bales e 


t al. 


1997: van Breueel et al. 


1998: 


Jarvis et al. 


2001 




De Breuck et all 2002: 


Willott et al. 


2003 


: Rocca-Volmerange et al 


. 2004: 


Brookes et al. 2006): HzRGs form a nar- 



row redshift sequence which traces the envelope of radio-quiet galaxies and is well modeled 
by the evolution of a stellar population formed at high redshift from an accumulation of up to 
10^^ Mq of pre-galactic baryonic material. The large-scale, double-lobed radio morphologies 
and enormous radio luminosities suggested early on that Hz RGs must have spinning super- 
mass ive black holes powering relativistic jets in their centers (|Reeslll978l : iBlandford &: Payne 
19821 ). With the more recent discovery that the stellar bulg e and central black hole masses 
of galaxies are closely correlated {e.g., iTremaine et al.ll2002l ). it is no longer a surprise that 
the parent galaxies of the most powerful radio sources occupy the upper end of the galaxy 
mass function. 

Despite two decades o f study since the initial discovery of the HzRG K — z rela- 
tion (jLilly fc Longairl Il984l ). the nature and tightness of the relation remains mysterious. 
The scatter in the relation is surprisingly low out to z ~ 1, and only increases mod- 
estly at higher redshifts as the observed fT-band samples rest-frame optical emission (e.g.. 
Fig. [H iDe Breuck et al.l |2002| ). Interpretations have generally relied on these observed 
near-IR observations probing the stellar populations of the distant radio galaxies {e.g., 
Rocca-Volmerange et al.ll2004j ). However, detailed interpretations have been complicated 
by (i) uncertain contributions of A GN-related light which are most important in the rest- 



frame UV {e.g., IVernet et al.ll200ll ). (ii) evolution of the host galaxy stellar population(s) 



and (iii) band-shifting, meaning that observed 2.2 /^m samples very different emitting wave- 
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lengths for different redshift sources, reaching into the rest-frame ultraviolet for the most 
distant HzRGs. 

In this paper, we pr esent observations o f a large sample of HzRGs obtained with the 
Spitzer Space Telescope (jWerner et al.ll2004l ). By observing the same rest-frame near- to 
mid-IR spectral range for sources over a large redshift range {1 < z < 5.2), we remove 
many of the complications which plagued previous studies. In particular, the 1.6 /xm peak 
of the stellar emission provides a reasonably robust measure of the stellar mass for stellar 
populations with ages ^ 1 Gyr (ISawickil l200ll ). While UV and optical emission has strong 
contributions from the youngest and hottest stars in a galaxy, near-IR emission primarily 
derives from the low-mass stars which dominate the stellar mass of a galaxy. Thus, while 
rest-frame UV and optical studies of galaxies are well-suited to probe galaxy star-formation 
rates, rest-frame near-IR studies are well-suited to probe galax y stellar rnasses . We discuss 
this approach and other recent results that might effect it (e.g. iMarastonI |2005| ) in detail in 
§6. 

Fig. [2] illustrates the advantage of rest-frame near-IR derivations of stellar masses. First, 
the impact of dust extinction falls sharply with wavelength, making quantities derived from 
rest-frame near-IR observations more than an order of magnitude less susceptible to uncertain 
dust extinction corrections relative to optical observations. Second, since the main sequence 
lifetimes of low-mass stars exceed the Hubble time, galaxy masses derived from rest-frame 
near-IR observations are relatively insensitive to the age of stellar populations or the star 
formation history. In particular, secondary bursts of star formation will affect UV and optical 
magnitudes significantly more than near-IR magnitudes. 

We here use mid-IR observations of a large sample of radio galaxies to probe their 
host galaxy stellar masses. By observing consistently on both sides of the 1.6 /zm peak 
of the stellar emission, we avoid complicated /c— correction effects and derive a reasonably 
robust measure of the stellar mass for stellar populations with ages ^ 1 Gyr. Observations 
at longer wavelengths allow us to determine the contribution of warm, AGN-heated dust 
emission to the rest-frame near-IR emission. Our paper is organized as follows. Section 2 
describes the Spitzer HzRG sample. Section 3 presents the Spitzer mid- to far-IR data 
and their reduction. Section 4 describes our fitting of the spectral energy distributions 
(SEDs) and §5 presents the results of this fitting. Section 6 describes how the derived 
rest-frame H— band luminosities are converted into stellar masses. Section 7 discusses and 
concludes this analysis. We present notes on individual sources as an Appendix. This paper 
presents our entire Spitzer data set and analyzes the bulk properties of the sample. Spitzer 
ob servations of individual sou rces from this program have b een the subject of detailed studies 
bv lVillar-Martm et all fl2006l . MRC 2104-242 at z = 2.49). lstern et al.N2006l . LBDS 53W091 
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at z = 1.55). lBroderick et al. I (120071 . PKS at z = 2.156) and lPe Breuck et all (120071 . 4C 23.56 
at z = 2.48). Throughout we assume a concordance models] of universe expansion, = 
1 — f^A = 0.3, fio = 1; and Hq = 70 kms~^Mpc~^ . Inferred luminosities presented in this 
paper are of the form uL^/Lq, where Lq = 3.9 x 10^^ W. 



2. The Spitzer High-Redshift Radio Galaxy Sample 



Our HzRG sample i s drawn from r adio galaxy surveys executed during the last 45 years 
BemiettI Il96lh. We haye sear ched both flux-limited surveys such 



starting with the 3CR: 



as th e 3C (ISpinrad et al 



20011 ). MG (IBennett et al 



98^ 



i yeCE dSales et al.lll997h . 7C Jhacv ejalJiigggl: 



1986: Lawrence et al. 



Wil 



1986[ ). and MRC jMcCarthv et al.lll996h. 



well as surveys filtered by their ultra-steep ra dio spectra (USS; e.g., [Chambers et al. 



ott et al. 



as 



1996 



Rottgering et al.l 119971 : iDe Breuck et al.l 120001 ) . The former provide a sample unbiased in 
their radio properties, but lack the high redshift {z > 2) sources mainly found in the USS 
samples. Because the USS samples also probe fainter radio flux densities, including such 
sources allows us to break the strong luminosity— redshift degeneracy of flux density-limited 
surveys. 

For the purposes of this work we define a HzRG as a radio galaxy above a redshift of 
one with a rest-frame 3 GHz luminosity greater than lO^^WHz^^. This choice is to ensure 
that we include only those objects with very powerful obscured AGNs. This value is similar 
to that used by other authors to separate AGN based up on their radio luminosities. The 



classical luminosity that separates iFanaroff fc Rileyl (Il974j ) type 1 and type 2 sources (which 
are morphologically distinct), converted from 178 MHz to 3 GHz, is 10^^'^ (10^^'^) WHz~^ 
assuming a steep (ultra-steep) radio spectral index of a = —0.8 (—1.3). A more recent 
radio- loud/radio-quiet division is Lsghz = 10^^-^ (lO^^-^)WHz 
upper limit to radio-quiet QSOs of LgcHz = lO^^ ^ (lO^^ i) WH 



(Miller et al. 


1990) and an 


was found bj 


I Greee et al. 



( 119961 ). both converted to 3 GHz luminosity assuming the same spectral index values. 



To examine how Spitzer-derived quantities depend on redshift and radio luminosity, we 
selected a representative subset of 69 targets from the parent sample of 263 known z > 1 
HzRGs (circa 2003; Fig. [3]). Table [1] lists the selected subset, which spans the full range of 



^These common values of the "concordance" cosmology underestimate luminosities (and hence stellar 
masses) by ^ 2 — 3% at z = 3 — 4 when compared to the latest Wilkinso n Microwave Aniso tropy Probe mea- 
surements, f^M = 1 — ~ 0.25, flo — 1, and Hq = 73 kms^^ Mpc~^ ( Spergel et al.l 12006 ). Lower-redshift 
HzRGs show an even smaller systematic adjustment, implying that the assumed cosmological parameters 
are inconsequential to our analysis. 
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redshifts fro m ^: = 1 to the redshift o f the most distant known radio galaxy, TN J0924— 2201 
at z = 5.19 fivan Breugel et al.lll999l ). Most sources were observed as part of Spitzer Cycle 1 
(Program ID number = 3329); a few observations derive from Guaranteed Time Observa- 
tions by Spitzer instrument teams. We required at least 15 sources in each Az = 1 bin, 
covering a range of rest-frame 3 GHz radio luminosity, Lsghz- To obtain a uniform set of 
radio luminosities of a sample of radio galaxies covering the entire sky, we obtained radio flux 
densitie s from the largest a ll-sky radio surveys, i.e., the 1.4 GHz NRAO VLA Sky Survey 
(NVSS: 



Condon et a 



19981), the 325/352 MHz Wester bork Northern Sky S urvey (WENSS: 



Rengelink et al.lll997l ) and the 365 MHz Texas Survey (jPouglas et al.lll996l ). These surveys 
all have ~ 1 arcmin spatial resolution, which has the advantage that the entire sample ap- 
pear essentially as point sources, thereby minimizing issues of missing flux from over-resolved 
emission. We chose 3 GHz as the fiducial rest frequency as it is the mid-point between the 
survey frequencies at the median redshift of our sample. This choice ensures that almost 
all radio luminosities are interpolated, rather than extrapolated as i s often required when 
deriving lower rest frequency luminosities {e.g., 151 or 178 MHz). iBlundell et al.l (119981 ) 
argue, however, that low frequency luminosities are a better measure of the total radio lu- 
minosity, as the emission at low frequency is dominated by the steep spectrum, isotropic 
radio lobes, while high-frequency emission may be contaminated by Doppler-boosted emis- 
sion from fiat-spectrum radio cores. Until the completion of all-sky <100 MHz surveys, such 
low frequency rest-frame luminosities cannot be uniformly derived for our sample. However, 
most HzRGs in our sample have weak radio cores, which a lthough flatter than the ultra-steep 



spectrum radio lobes, still have steep radio spectra (e.g.. lAthreya et al.lll997l ). and thus are 



not expected to significantly boost the rest-frame 3 GHz luminosity. 

We have preferentially selected sources with the most supporting data. Forty-seven 
of the 69 HzRGs in our sample have HST data, representing over 1.5 Ms of HST time. 
We also preferentially included sources previously selected for guaranteed time observations 
(GTO) by the Spitzer instrument teams, thereby minimizing our Sj^zteer observatory request 
while maintaining a uniformly selected sample. Finally, we preferentially selected sources 
with millimeter and sub-millimeter observations in order to have information on the SEDs 
longward of the Spitzer coverage. Note that the sample does not suffer a bias toward cold 
dust properties as we select only on the basis of (sub-)mm observations, not detections. 



3. Infrared Observations and Reductions 



The Spitzer observations consist of imaging photometry in eight bands, using the entire 
instrument complement of the observatory. Data from the Infrared Array Camera (IRAC; 
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Fazio et al.ll2004l ) comes from all four bands (3.6, 4.5, 5.8 and 8 /im — synonymously referred 
to as channels 1 to 4). The entire sample of 69 HzRG s was observed with IRAcj^. Data from 
the I nfraRed Spectrogra ph (IRS; iHouck et al.l 120041 ) obtained in "peak-up imaging" mode 
(c./., iTeplitz et al.ll2005l ) provides flux densities at 16 /im for the 46 objects in our sample 
at z > 2. For the lower redshift portion of our sample, IRAC 8 /im images adequately cover 
the SED longward of the 1.6 p,m stell ar bump. Finally, the Multi-band Imaging Photometer 
for Spitzer (MIPS; iRieke et al.l |2004| ) was used to image a subset of our sample using all 
three available bands (24, 70 and 160 /xm). Due to the uncertain ability of MIPS to image 
against the Galactic infrared background at the time of our GO Cycle 1 submission, only 
the 26 HzRGs from our sample lying in "low" background regions (5*24^111 < 20MJysr~^) 
were observed with MIPS. Tabled] presents the 69 HzRGs comprising our sample, along with 
their coordinates, redshifts, rest-frame 3 GHz radio luminosities, 64kpc i^'-band magnitudes 
(Vega) and which Spitzer instruments they have been observed with. 

The 141 unique target/instrument observations comprising our survey were released 
from the Spitzer archive between Octob er 2004 and June 2006. Fourteen observations were 
part of the initial GTO programs {e.g., IStern et al.ll2006l ). In general, our data reduction 
begins with the basic calibrated data (BCD) output of the Spitzer Science Center pipeline, 
using the most current version number at the time the data were initially released [e.g., 
S9 to S13). Improvements over the generations of pipelines implemented here were minor 
and did not warrant reprocessing data. The most useful augmentation was the addition of 
a global world coordinate system to the IRS peak- up images in Sll. Note that the units 
of Spitzer BCD images are in surface brightness units of MJy ster~^, where 1 MJy ster~^ 
corresponds to 23.5 /xJy per square arcsecond. 



3.1. IRAC Observations 

Most of the IRAC observations consisted of four dithered 30 s exposures in each of the 
four IRAC channels. Fo r each channel, the BCDs were mosaiced using the MOPEX package 



(IMakovoz fc Khanll2005l ) from the Spitzer Science Center and re-sampled by a factor of two 
to give a pixel scale of 0''61. The MOPEX outlier {e.g., cosmic ray, bad pixel) rejection was 
optimized for the regions of deepest coverage in the center of the maps, corresponding to 
the location of the targeted HzRGs. The eight GTO IRAC observations were all deeper 
than our standardized program and were reduced in a similar manner with MOPEX. Source 



^Although 3C65 was observed in only two bands due to an astrometric error in NED which we have since 
had corrected. 
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extraction was performed with SExtractor ( iBertin &: Arnoutall996l ) in dual image mode 



using the 3.6 fim channel for source detection. All IRAC BCD images have astrometry 
derived from 2MASS survey and are good to < 0'.'5 so we can be confident of the radio 
galaxy identifications. 

We generally used the a 7" diameter aperture for photometry, though in a few cases 
where the source was faint and/or had a close companion, we extracted photometry using a 
3''5 diameter aperture. The aperture corrections, applied after the conversion from surface 
brightne ss to flux density , are listed in Table [2] and are taken from the Spitzer First Look 



Survey (iLacy et al.l l2005l ). The resulting flux densities are listed in Table [3] and include 
the total uncertainty in the IRAC photometry; i.e., the statistical uncertainties are added 
in quadrature to the systematic uncertainty of ~ 10%. The systematic uncertainty is the 
total uncertainty from several different effects: the uncertainty in the IRAC flux density 
calibration, a fiat fielding color dependent uncertainty and variation of the PSF across the 
field of view. Non-detections have 3cr upper limits listed. The statistical uncertainties in 
the mosaiced images were calculated from the pixel-to-pixel rms from single BCDs (with 
the native pixel scale), converted to /xJy, multiplied by the square root of the area of the 
aperture in square pixels, divided by the square root of the number of frames, and multiplied 
by the appropriate aperture correction. In all cases, the radio galaxies are well detected in 
the 3.6 and 4.5 /im channel^. In most cases, the radio galaxies are also detected in the 5.8 
and 8 fim channels. In no cases were the HzRGs resolved in channel 1, the highest resolution 
of all the Spitzer imaging bands with a PSF FWHM of 1''44. 



3.2. IRS Observations 

The IRS observations consisted of just two nodded exposures of the 16 /xm peak-up 
imager, each of 61s duration. The observations were executed in Spitzer Cycle 1, before 
IRS "peak-up imaging" was offered as a supported observing mode. However all the data 
were reduced with the more recent S13.2.0 version of the archive pipeline, which provides 
astrometry to the BCD products accurate to 0''5 relative to the 2MASS reference frame. 
Sub-images were made of just the blue (16 /im) "peak-up" field of view in which the radio 
galaxy was centered. Difference images were then created by subtracting one frame from the 
other, thereby removing background emission and latent patterns in the array. 



■^3065 is not observed in the 4.5 /im channel. 6C0140+326 is probably detected at 3.6 and 4.5 //m, but 
is blended by a foreground galaxy. Likewise 3C294 is detected in all four channels, but the proximity of a 
bright star means only upper limits can be determined for channels at 3.6 and 4.5 /xm 
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The dithered difference images were then combined with bad pixels and ffaring pixels 
masked. Flux densities were extracted with the IRAF DAOPHOT task using the parameters 
listed in Table [H Random uncertainties were determined from the rms in each BCD and 
converted to /xJy. The total uncertainties in this table include the random and the 7.8% (6% 
zero point uncertainty and 5% color correction) systematic uncertainties. The systematic 
and random uncertainties are combined in quadrature and used for the SED fitting discussed 
in §4. These flux densities and uncertainties are listed in Table [31 Non-detections are 
represented by 3a upper limits. All 46 of the Spitzer HzRG sample at z > 2 were observed 
with IRS, as was MRC 0211-256 {z = 1.300) due to an error in the submitted AOrA We 
have robust 16 /xm detections of the radio galaxy in 35 of the 47 observation^ 



3.3. MIPS Observations 

The MIPS observations were restricted to the 26 HzRGs with "low" mid-IR background 
< 20MJysr^^) and consisted of one cycle of 30 s exposures at 24 /xm, three cycles of 
10 s exposures at 70 /im, and five cycles of 10 s exposures at 160 fim. The 24 fim BCDs were 
median filtered to remove a regular stripe artifact pattern prior to mosaicing. Both the 24 
and 70 fim data were then mosaiced with MOPEX using standard inputs and with the pixels 
re-sampled by a factor of two. Twenty-two of the 24 fim observations resulted in detections, 
but only four radio galaxies are detected at 70 /im. None of our HzRG sample are detected 
at 160 ^m. As we shall see below, the upper limits at long wavelengths are still quite useful 
for constraining the mid-IR SED. Flux densities were extracted with the IRAF DAOPHOT task 
using the parameters listed in Table [2l For the 70 /im band we use an aperture correction 
based on a steep spectrum {Fx oc A°). Again, Table [3] provides flux densities and upper 
limits (3(j) for all three MIPS cameras. 

The random uncertainties (and 3cr upper limits in Table [3]) are derived from the BCDs 
and converted to yuJy. The color correction can be quite significant for these broad MIPS 
bands which were calibrated using Galactic stars, e.g., typically for a Raleigh- Jeans spec- 
trum, F\ oc A~^. For the 70 and 160 /im bands we include a 5% systematic uncertainty to 
account for the uncertain color correction. For 24 /xm we typically find that the SED is very 
non-stellar. As a first approximation, we apply the color correction by multiplying by 0.961 



''It was inadvertently confused with USS 0211-122 (z = 2.336). 

^Except for the non-detection of the highest redshift radio galaxy, TN J0924— 2201 (z = 5.195), the 16 /ini 
non-detections are evenly distributed across z 2 — 4. 
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(from the MIPS data handboolo) for XFx = constant, which approximates the observed 
SEDs of our HzRGs. We further include a 2% systematic uncertainty for the remaining 
unknown cofor correction. The final systematic uncertainties are are 4.5%, 8.6% and 13% 
for 24, 70 and 160 /xm, respectively, including flux calibration uncertainties from the Spitzer 
Observers Manual and color correction uncertainties described above. The systematic un- 
certainties are combined in quadrature with the statistical uncertainties and are presented 
in Tabled 



3.4. Mid-IR Colors of HzRGs 



Recent work has shown that IRAC mi d-IR colors reliably isolate spectr oscopically 



confirmed, luminous, unobscured AGN (e.g.. iLacy et al.ll2004l : IStern et al.ll2005l ). However 



since the surface density of sources with mid-IR colors similar to unobscured AGN is much 
higher than that of optically-selected (type 1) quasars, it has been inferred that obscured, 
type 2 quasars are likewise identifiable from their mid-IR colors. Due to their optical faint- 
ness, type 2 quasars are not as well studied as their bright, type 1 cousins. However, they 



are a natural consequence of AGN unification schemes (e.g.. lAntonuccilll993l ). and models of 



the h ard X-ray backgroun d imply that they outnumber type 1 quasars by factors of several 



e.g. 



Treister et al.l 120041 ). 



Our Spitzer HzRG sample represents the largest sample of confirmed type 2 quasars 
with IRAC photometry obtained to date. Fig. H] presents the measured mid-IR colors of the 
49 HzRGs with detections in all four IRAC bands and the nine HzRGs with two-band IRAC 



detections on the color-color plots of iLacy et al.l (120041 ) and lStern et al.l (120051 ). As expected, 
the HzRG samp l e tend s to fall within the AGN wedges identified in each of these diagrams. 
The iLacy et al.l (120041 ) re gion includ e s nea rly the entire HzRG sample: 51/58 (88%) of 



HzRGs plotte d match the 



the criteria of iLacy et al 



Lacy et al.l (120061 ) c olor criteria ( which are slightly updated from 



20041 ) . The smaller IStern et al.l (120051) region misses a few more 



HzRG, with only 41/58 HzRGs (71%) matching the IStern et al.l (120051 ) color criteria. For 
both plots, nearly every HzRG is within la of matching the AGN color criteria. 

It should be noted that while the obscured and unobscured AGN generally lie within the 
"AGN wedges" in Fig. Hj, these two classes of luminous AGN are presumed to have different 
energetics powering their observed mid-IR, IRAC-band SEDs. Luminous, unobscured AGN 
tend to have power- law SEDs due to the central engine across the mid-IR regime [e.g., 



Richards et al 



20061 ). This leads to rising spectra through the IRAC passbands and places 



' jhttp : //ssc . spitzer . caltech.edu/mips/dli/inipsdateLheLnd.book3 . 2.pd.f 
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such sources squarely in the wedges at most redshifts. The caveat is that strong features 
can cause excursions as they redshift into the passbands; e.g. Ha at 2; 4 is in the 3.6 fim 
band (e.g., Bian et al., in prep.). In contrast to this, we presume that starhght plays a 
dominant role in the shorter wavelength IRAC passbands for most obscured AGN, with 
AGN-heated dust dominating the reddest IRAC passbands. Indeed, a clear dip between 
the stellar component and the hot dust component is seen in many sources in our sample. 
Consider, as an example, the left-hand panel of Fig. 3. The redshifted galaxy tracks have 
[3.6] — [4.5]> 0.2 (Vega) at z > 1 due to the 1.6 /im peak of the stellar SED red-shifting into 
the 4.5 yum IRAC passband. The red [3.6] — [4.5] colors is then usually due to the z > 1 stellar 
populations while the red [5.8] — [8.0] colors are due to the AGN. The next section will concern 
using the longer- wavelength Spitzer data to model and subtract the AGN contribution to the 
mid-IR SEDs of our sample, and thus to derive the rest-frame near-IR stellar luminosities. 



4. Modeling of the Rest-Frame Near-IR and Mid-IR 

We begin our analysis of the broad-band IR SEDs by considering the galaxies with the 
most extensive Spitzer coverage. This smaller sample of 26 galaxies illustrates the shape and 
range of HzRG SEDs in the IR (§4.2), and is then used to inform analysis of the sources 
with more restricted Spitzer imaging (§4.3). 

Despite the extensive near-IR data available for many of the HzRGs, we elect not to 
include these data in the SED fitting (though the i^-band photometry is over-plotted in 
the SED plots). This omission has several reasons and, as described in §5, should have 
minimal impact on the derived stellar masses, the primary goal of our analysis. Foremost, 
we omit the near-IR data as it samples the rest-frame UV to optical wavelengths of our high- 
redshift galaxies. At these wavelengths, several issues arise. First, scattered and direct AGN 
continuum is blue i n color and will thu s be a larger contributor to the SED at wavelengths 



below 4000 A (e.g.. IVernet et al.ll200ll ). Second, high equivalent widt h emission lines can 



signifi cantly affect the near-IR photometry for certain redshifts [e.g., IVillar-Martm et al. 



20061 ) . Indeed, when Ha is in the bluest IRAC channels (for our highest redshift sources), 
we do not include the affected channel in the SED modeling. Third, possible contributions 
from younger stellar populations will be most pronounced at shorter wavelengths. Fourth, 
rest-frame optical and UV data are more susceptible to dust extinction than rest-frame 
near-IR data. The final reason to omit the near-IR data is the large inhomogeneity of 
the data quality, depth, and analysis techniques. We are generally relying upon published 
photometry done on 3 m to 10 m telescopes over the past two decades and extracted with a 
range of aperture sizes. Only in a few cases are we able to obtain the images from telescope 
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archives. 



For the sake of homogeneity of the analysis, the SED modehng described below relies 
solely upon Spitzer photometry. In a forthcoming paper, we study the detailed stellar popula- 
tion properties of the HzRG sample, combining Spitzer data with archival data and data from 
our ongoing observing campaigns at Palomar Observatory and the VLT (A. Rettura et al., 
in prep.). For a few galaxies, we have already performed the full analysis using all available 



(opti cal to mid-IR) data {e.g., IVillar-Martfn et al.ll2006l : IStern et al.ll2006l : iDe Breuck et al. 
20071 ). As shown in §5, the derived stellar masses are consistent with the results of the 
simpler analysis presented here. 



4.1. The Model 

Twenty-six of our radio galaxies have both IRAC and MIPS observations, of which 
only four are not detected at 24 /zm (or in the longer wavelength MIPS channels). These 
four sources include three of the highest redshift radio galaxies and one radio galaxy at 
z ~ 1.5. In the latter case, this non-detection may be indicative of silicon absorption at 
rest-frame 9.8/xm. Alternatively, t he undetect e d z ry 1.5 radio galaxy, LBDS 53W091, is 



a low luminosity radio source and IStern et al.l (120061 ) argue that its mid-IR non-detection 



likely indicate a SED dominated by old stellar populations and not by hot, AGN heated 
dust. One source, 6C 0140-1-326, cannot be distinguished from a bright foreground galaxy 
at 3.6 and 4.5 fim so we do not use it in our modeling here. Hence our final sample numbers 
21 HzRGs with IRAC and MIPS detections. 

For the 21 HzRGs with MIPS detections, we fit the SEDs with toy models comprised of 
a Composite Stellar Population (CSP) and multiple black-body (BB) dust components. Our 
primary goal is to simply derive the rest-frame near-IR stellar luminosity, which can then 
be translated into a stellar mass assuming reasonable stellar population synthesis models 
(§5). We are also, in these cases with MIPS detections, able to estimate the rest- frame 5 /xm 
monochromatic luminosities, a possible measure of the obscured AGN power (jOgle et al. 



20061 ). The difficulty is that in many sources a steeply rising continuum is visible in the 
redder IRAC bands, indicative of the tail of hot, AGN- heated dust contaminating the rest- 
frame near-IR emission. Owing to the sparse sampling of the SEDs at these wavelengths, our 
goal is not to fully explore the astrophysics of the continuum emission at these wavelengths, 
but to use the MIPS (and IRS, where available) photometry to model the dust SED and 
subtract its contribution at rest-frame near-IR wavelengths. We do not apply Galactic 
extinction corrections as they are negligible (< 0.5%) at the Sj^ifeer wavelengths. Also, since 
emission lines longward of Ha are not strong enough to significantly affect our modeling 
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e.g., iHaas et al.ll2005l : lOgle et al.ll2006l : ICleary et al.ll2007f ). they are ignored in the analysis. 



Additional systematic issues in our modeling are discussed in §4.4 and §6.3, but we stress 
here that the stellar masses are very insensitive to the choice of BBs for the longer wavelength 
component. 

We model the radio galaxy IR SEDs with a four component model consisting of an 
elliptical galaxy CSP template and three BBs of different dust temperatures. Best-fit pa- 
rameters for the models are determined using standard minimization techniques. The 
uncertainties on the flux densities used to determine the of a specific fit include both ran- 
dom (from the observations) and systematic {e.g., the detector zero points and the shape 
of the SED) uncertainties, which are added in quadrature (see §3). The full fitting is only 
done for sources with MIPS 24yum detections as we need both enough data points to fit our 
model and a detection longward of the IRAC bands to provide sufficient leverage on the dust 
component (see §4.2). For sources without MIPS detections or observations we do not have 
enough data points to perform a secure fit of this four component model. Hence in many 
cases we only get an upper limit to the stellar component, though in some cases we can get 
a likely value for the stellar luminosity when the IRAC bands fit the stellar component well 
(see §4.3). 



4.1.1. Stellar Component 



The CSP t emplate used in our SED model ing consists of an elliptical galaxy modeled 
with PEGASE . 2 (Fioc fc Rocca-Volmerangdll997l ). The adopted star formation history (SFH) 
scenario has been shown to result in reliable photometric redshifts with the code Z-PEG 
( iLe Borgne fc Rocca-Volmerangell2002l ) and to well reproduce early - type galaxy SEDs (from 
rest-frame UV to near-IR) up to 2; ~ 1.25 (IRettura et al.l l2006l ). Specifically, the star- 
formation rate is proportional to the gas density, SFR cx pgas, and the astration (recycling) 
rate is equal to 3.33 Gyr^^. Gaseous exchanges with the interstellar medium are represented 
by infall and galactic winds. The former process simulates mass growth and the latter 
represses any star-formation after the first Gyr by dissipating the gas from the galaxy and 
hence driving the gas fraction towards zero. The reader is referred to the aforementioned 
papers for m ore details on the adopted templ ate parameters and the star- formation histories. 
Additionally iRocca-Volmerange et al.l (120041 ) demonstrates how this elliptical galaxy model 
works for the most massive g alaxies that po pulate the lower envelope of the K — z relation. 
For the current work, we use a lKroupal (1200 ll ) initial mass function (IMF) and dust-free model 
templates. We note that in the templates used here, the metallicity is not a free parameter; 
it evolves consistently with the SFH, reaching solar values at the ages of interest in this 
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study. We assume a galaxy formation redshift of Zform = 10, and we discuss the systematic 
effects of this choice, along with the other details of the galaxy models in §6.3. In brief, 
however, though the derived stellar masses vary somewhat with the model choices (mainly 
by a simple offset depending on the choice of IMF), we find that the derived luminosities are 
relatively robust against model uncertainties for stellar population ages greater than 1 Gyr. 



4- 1-2. Dust Component 



The longer wavelength modeling assumes three dust components, each of a different 
temperature. The coldest of the BBs is chosen to have a temperature 60 K, representativ e 
of the large cool dust reservoirs associated with radio galaxies {e.g., iGreve et al.l l2006l ). 
although this component does not contribute to the rest-frame near-IR. The second BB 
component has a fixed temperature of 250 K. For the final, hotter BB component we allowed 
the temperature to vary from 500 to 1500 K (in steps of 50 K) and interpret it as being 
associated with hot dust heated by the central AGN. This upper limit of 1500 K was chosen 
as dust sublimates at this temperature. Systematic effects from the chosen BB temperatures 
are discussed in §6.3; in brief, we find they have little effe ct on the derived stellar masses. We 
have also experimented with toroidal AGN models (e .g., iPier &: Kroliklll992l : iGranato et al. 



19971 : iNenkova et al.ll2002l : ISiebenmorgen et al.l 120041 ). However, the paucity of our infrared 
data, consisting simply of photometric points and not infrared spectra, makes fitting with 
detailed models ill-constrained. We emphasize that this dust modeling is not completely 
physical; however, the modeling is sufficient for the primary goal of this paper, which is 
to simply estimate and subtract any contribution by (hot) dust to the rest-frame near-IR 



emission. 



4.2. Fits to HZRGs with MIPS Detections 

The best fit models for the 21 HzRGS with MIPS detections are shown in Fig. [51 Ob- 
served wavelength is indicated along the bottom axis and rest wavelength is on the top axis. 
Open circles, not used in the fitting, are either (i) K-hand data or (ii) potentially contami- 
nated by strong emission lines {e.g., Ha). Upper limits are indicated by downward arrows. 
We separate the contribution of host dust and stellar emission at rest-frame H— band, which 
can then be used to derive stellar masses (§5). Table H] presents the total if— band luminosi- 
ties, the stellar H— band luminosities, the stellar fractions at H— band, stellar masses, and 
rest-frame 5 /xm monochromatic luminosities for sources with full SED fitting. This fiducial 
mid-IR wavelength, 5 /xm, was selected as it is always short-ward of the 24 /im MIPS obser- 
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vations, even for our most distant HzRGs with 24/xm detections [z ~ 3.8). The formal 
uncertainties (95.4% confidence hmits) of the stellar luminosities from the fitting are listed 
in Table m They are generally of the order ~ 25% and are mainly due to the uncertainty in 
the IRAC flux densities, which, in turn, are dominated by systematics. 

4.3. Fits to HzRGs without MIPS Observations or Detections 

We can constrain, or in some cases quantify, the stellar luminosities for sources without 
MIPS observations or detections using the results of the full SED fitting described above. 
For many HzRGs we clearly observe the if— band stellar peak (see Fig. Even in those 
cases where we do not see the if— band peak, IRAC channels one and two are typically 
dominated by the stellar component. At the least, we can obtain an upper limit to the 
stellar mass by demanding that the stellar component does not exceed any of the IRAC 
photometric points. Hence we can obtain upper limits to the stellar masses of these HzRGs. 
Typically, this upper limit is restricted by just IRAC channel 1, as the observed mid-IR 
SEDs are often quite red. In many cases, the maximum fit is constrained by IRAC channels 
1 and 2, and occasionally by longer wavelength channels. In almost half of the radio galaxy 
sample, the stellar SED is well fit in two or three channels (those at shorter wavelength), 
so we nominally fit the SED and derive stellar luminosities. The results of the fitting for 
the HzRGs without MIPS observations or detections are shown in Fig. El which adopts the 
same symbol convention as Fig. [5] bar the inclusion of BB dust components. As can be seen, 
the 16 /im observations, where present, generally show a strong excess above the maximally 
fitted stellar template. In what follows, we distinguish those HzRGs whose IRAC SEDs are 
well fit by a stellar population (e.g., where multiple IRAC channels are well fit by the stellar 
SED model) from those where the IRAC photometry has merely provided an upper limit to 
the rest-frame near-IR stellar luminosity. 

4.4. Systematic Effects in the SED Modeling 

We discuss here two possible systematic effects in the SED modeling: the adopted BB 
temperatures and potential photometric offsets in the longer wavelength photometry due to 
strong spectral features redshifting into the Spitzer passbands. The latter effect mimics the 
effects of alternate SEDs with strong absorption and/or emission features. As seen below, 
neither systematic issue strongly affects the derived rest-frame near- and mid-IR luminosities. 

As the cooler dust components have negligible influence at rest-frame ii— band, the 
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exact temperatures used are not vital to the final derived stellar luminosities and masses. 
Modeling observations of 3C radio gala xies obtained with the Infrared Space Observatory, 
Rocca-Volmerange fc Remazeilled (120051 ) find the far -IR SEDs are wel l fit w ith two BB com- 
ponents of temperatures 40 ± 16K and 340 ± 50 K. [Archibald et al.l (120011 ) . modeling sub- 
millimeter observations of HzRGs, likewise adopts an optically-thin, isothermal grey-body 
of temperature 40 K. When we redo our analysis adopting this cooler cold dust compo- 
nent rather than the previously adopted 60 K, the derived rest-frame near/mid-IR lumi- 
nosities change by less than 0.1%. Li kewise, using the 340 K warm dust component of 
Rocca-Volmerange fc Remazeilled (120051 ) instead of the previously adopted 250 K, we find 
that the derived rest-frame near-IR luminosities change by less than 0.1%, although the 
mid-IR luminosities increase by ~ 40%. 

The fitting algorithm we adopt assumes a very simplistic m odel of the long- wavelength 
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even with a line width of 1. 4 ^m and r = 1.5 we would only observe a 30 — 40% decrease 



in the observed flux density. I Ogle et al.l (120061 ) also note that there are indications that the 
power-law continuum is truncated at the sublimation temperature. 

We have investigated the effects of these potential features by re-running the fitting 
procedure, varying the 16/im and 24/im flux densities to simulate those passbands falling on 
strong emission or absorption features. We varied the 16/im flux densities by 10% to account 
for weaker PAH features at 6.4/im passing through the IRS bandpass. We varied the 24/im 
flux densities by 50% to account for both the strong 7.7/im PAH emission feature and the 
strong 9.8/im silicate absorption feature. These experiments only changed the mean 5 fim 
luminosities by ~ 5% (for varying the 16 fim flux density) and ~ 25% (for varying the 24 fim 
flux density). Given the lack of observational evidence for these features so far we consider 
these estimates conservative. The effects on the mean stellar luminosity (and hence mass) 
are ^ 2% in both cases. 



5. Results of SED fitting 

5.1. Stellar Fractions 

Figures [7] and [H] show the fraction of the HzRG rest-frame H— band luminosity which 
is modeled as stellar in origin, plotted against total H— band luminosity, stellar H— band 
luminosity, redshift, radio power, and rest-frame 5 fim luminosity The stellar fraction is 
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defined as the ratio of the flux densities found from convolving the stellar and total SEDs 
with a H— band filter curve. 

Only the 21 HzRGs with MIPS detections are presented in these figures and their 
iJ— band luminosities (stellar and total), stellar fractions, inferred stellar masses (see §6), 
and monochromatic 5 /im luminosities (see §5.4) in Table HI As can be seen, half the sample 
has the rest-frame near-IR luminosity strongly dominated by starlight, with > 95% of the 
emission apparently stellar in origin. Another 20% of the sample have stellar emission ac- 
counting for 60 — 90% of the rest-frame near-IR flux, and thus still dominating the emission 
at this wavelength. There appears to be a weak inverse correlation between stellar fraction 
and total H— band luminosity whilst the stellar fraction does not vary with stellar H— band 
luminosity. This suggests that the underlying host galaxy population is homogeneous. No 
evolutionary trends in the stellar fraction are evident, nor does the stellar fraction appear to 
correlate with either radio power or rest-frame mid-IR luminosity. This result is expected for 
the following reasons. The radio emission at rest-frame 3 GHz is dominated by isotropic lobe 
emission. Likewise, rest- frame mid-IR samples largely isotropic emission coming from near 



the central engine {e.g., iHeckman et al.lll994j ): at these wavelengths, the emission should 



be largely immune to extinction. The stellar fraction at H— band, however, is very depen- 
dent upon obscuration of the AGN wh ich, in the currently preferred unified schemes {e.g., 



Antonuccilll993l : lUrry &: Padovanilll995l ). is a function of orientation. AGN with little obscu- 
ration of the central engine will have higher rest-frame near-IR luminosities, show broad lines, 
and have low stellar fractions. Indeed, the three sources in our MIPS sample with the highest 
rest-frame H-hand luminosities, PKS 1138-262, TX J1908+7220 and B3 J2330+3927, both 
show AGN signatures in the near-IR (see Appendix) and are indicated as upper limits in Ta- 
ble H] and the relevant figures (Figs. [TJ [8]). The stellar fraction is expected to be an indicator 
of orientation and should thus not correlate with the orientation-independent parameters. 
In a future publication (Seymour et al. in prep.), we shall analyze how the stellar fraction 
correlates with properties which are expected to depend on orientation, such as radio core 
fraction and radio spatial extent. 



5.2. Rest-frame Near-IR Stellar Luminosities 

The rest-frame H— band stellar luminosities are found mainly to lie in the 10^^ — IO^^Lq 
range, as illustrated in Fig. O Overlaid are PEGASE models illustrating the luminosity of 
an elliptical galaxy model with Zform = 10 for masses of 10^^ and IQ^^Mq. The derived 
HzRG stellar luminosities are consistent with extremely massive hosts, with masses of a 
few xIO^^Mq out to z = 4. Derivations of masses and their associated uncertainties are 
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presented in §6. 



This result confirms earlier work on the masses of radio galaxies. iRocca-Volmerange et al 



( 120041 ) find that most radio galaxies have K— band magnitudes consistent with stellar masses 
of 10^^ — lO^^M© with a few having magnitudes implying masses > lO^^M©. These sources 
with bright K— band magnitudes may have emission lines or scattered AGN light contribut- 
ing to them so their K— band should be viewed as upper limits. 

We also plot the rest if— band stell ar luminositi e s of a sample of sub-mm galaxies 
(SMGs) in the GOODS North region from lBorys et al.l (120051 ) as red crosses in Fig. [91 The 
H — band luminosities are re-derived using our modeling technique wit h the data pr e sente d 



m 



Borys et al.l (120051 ). Our analysis produces similar luminosities to iBorys et al.l (120051 ). 



However, due to different assumed star-formation histories and IMFs, our derived stellar 
masses systematically differ by 25%. These SMGs have a similar range of rest-frame H— band 
luminosities to the HzRGs, implying similar stellar masses. There is a probable evolutionary 
link between SMGs and HzRGs, though the number of similarities between these populations 
is comparable to the number of differences. On the one hand, SMGs have a much higher 
areal density than HzRGs: the former are relatively numerous in small-area (100 arcmin^) 



surveys [e.g. 



Borys et al.ll2003l ). while the latter are much rarer and typically herald from 
large-area or full-sky radio surveys. On the other hand, man y radio galaxies are strong 



sub-mm emitters {e.g., [Archibald et al.ll200ll : iGreve et al.ll2006l ) and would qualify as SMGs 



irrespective of their strong radio emission. 



5.3. Mid-IR Luminosities 

The rest-frame 5 /xm monochromatic luminosities presented in Table H] are generally all 
above 10^^ Lq. In Fig.[TO]we show the distribution of the 5 /im luminosities with redshift. The 
solid symbols and arrows indicate the 21 HzRGs with MIPS detections and the five HzRGs 
with MIPS upper limits. The open symbols indicate the 13 HzRGs with IRS observations 
but no MIPS observations where the IRS band is close to rest-frame 5 /im. These IRS-derived 
luminosities are calculated directly from the 16 /xm flux densities, and thus do not account 
for the shape of the SED or the related offset between rest-frame 5/im and the rest-frame 
wavelength probed by the 16/im observations. The uncertainty on the IRS-derived rest-frame 
5/im luminosities is estimated to be < 50%. 

The high-luminosities, Ls^m > 10^^ L©, in Fig. [TO] imply that the HzRGs are re- 
radiating a substantial amount of energy, most likely re-processed emission from the ob- 
scured AGN heating the surrounding dust. The imphed total (5 — 1000/im) IR luminosi- 
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ties will be greater than the mid-IR luminosities presented here and hence most of these 
HzRGs would be classified as LIRGs (-^^s-ioooyum > lO^^i^o), or in many cases as ULIRGs 
(L5_iooo/im > 10^^ Lq) had they been initially identified with no kno wledge of their ra - 
dio properties. These luminosities cover a similar range as found by lOgle et al.l (120061 ) 
for a sample of quasars and radio galaxies from the 3CRR catalogue: z/Ly(rest 15 /im) 
= 10^'*"*^^ erg s~^. Luminosities this high confirm that powerful radio galaxies harbor "hidden 
quasars", i.e., super-massive black holes accreting ne ar the Eddington limit. A few sources 
have lower luminosities, classified by lOgle et al.l (120061 ) as being in the mid-IR "weak" regime, 
z/L;y(rest 15 yum) < 10'^^ ergs"^, approximately equivalent to z^Li^(rest 5 /im) < 10^^ Lq, in- 
dicating that they probably have lower accretion rates, and are possibly non-thermal, jet- 
dominated sources. 



The alternate explanation for these high mid-IR luminosities, star-form ation, is unlikely. 
Desp i te the high inferred s tar-formation rates from the sub-mm observations (lArchibald et al. 
200ll : iReuland et aPboOsh . > 500MqYt-^ for the ~ 30% detected and < 5OOM0yr-i for the 
non-detections, the star-formation is widely scattere d as seen in the rest-frame UV/qptical 



images which show widely distribu ted morphologies (jMcCarthy et al.lll987l : iBest et al.lll996 



Zirm et al.ll2003l : iMiley et al.ll2006l ). Hence the star-formation is unlikely to be dense enough 
to reach temperatures where it would dominate at rest-frame 5 /xm. 



5.4. Near- to Mid-IR Luminosity Ratio 

Figure [TT] presents the stellar near-IR to total mid-IR (-Z^H*^V-^5°/im) rest-frame luminos- 
ity ratio plotted against both redshift and 3 GHz luminosity. The ratios range from ~ 0.1 
to ~ 6 and we find that the near- to mid-IR luminosity ratio gradually decreases with both 
increasing redshift and with increasing radio luminosity. If, as discussed earlier, the mid- 
IR luminosities are indicators of the power of the AGN [i.e., the AGN accretion rate at 
the epoch of observation), then for a given host stellar mass, the accretion rate apparently 
declines at lower redshift and lower radio power. This correlation suggests evolution with 
redshift, but may also be due to selection: our mid-IR luminosities are based on the 24 /xm 
detections which are fewer at higher redshift. The trend with radio power is perhaps more 
interesting, showing that HzRGs with more luminous radio emission are also accreting at a 
higher rates (for a given stellar mass) . This result is also c onsistent with the radio luminosity 



emission line luminosity correlation (jWillott et al.lll999l ). This trend with redshift is prob 



ably also related to the slight trend, albeit with a large scatter, when directly comparing the 
mid-IR to radio luminosity (see Fig. [12]), although that trend is possibly due to the scahng 
of luminosities with distance. 
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Photometric Stellar Masses 



The stellar luminosities shown in Fig. [9] are converted to masses using the mass-to- 
light ratio^ derived from the PEGASE . 2 early-type galaxy model with ^form = 10 described 
in §4.1 and are presented in Fig. [131 The plotted errors are from the SED fitting of the 
luminosity; systematic errors to the stellar masses are discussed in §6.3. Hence we show that 
the stellar masses of the hosts of radio galaxies generally have masses distributed between 
10^^ — 1O^^'^M0 with remarkable little scatter across the redshift range of 1 < z < 4. The 
one HzRG with larger uncertainty a.t z = 5.2 suggests this trend may hold to higher redshifts 
still with the implication that this source had to have formed early and quickly. Given the 
lifetimes of radio galaxies are short, 10 — 100 Myrs, these sources are likely to be representative 
of a larger population of massive galaxies at these redshifts. Massive g^alaxies are now being 
found in great numbers across 1 < -2 < 
van Dokkum et al.ll2006l ). 



3 {e.g., iLabbe et al.l l2005l : iPapovich et al.l 12006 



We compare our HzRG sample to SMGs studied by iBorys et al.l (120051 ) (red crosses 
in fig. [T3|) . HzRGs and SMGs are amongst the most massive objects over this redshift 
range, showing that both classes of galaxy probe the most massive g alaxies in th e high- 
redshift universe, though we note that HzRGs are much rarer beasts. iBorys et al. I (I2OO5) 



studies 13 SMGs identified in the GOODS survey. For comparison, the most luminous radio 
source known in the GOODS survey i s VIA J123642+621331, a 470/iJy source at z = 4.24 
identified by (IWaddington et al.lll999l ). With Lsghz ~ 10^^ WHz~\ the radio emission from 
this galaxy is most likely comi ng from an AGN, no t star formation, as suggested by European 
VLB! Network observations (iGarrett et al.ll200ll ) which show the source to be unresolved. 
Nevertheless, it's radio power does not qualify it as a HzRG according to our selection criteria 
defined in 52. 



6.1. Radio Power Dependence on Stellar Mass 

Fig. [H] presents the derived stellar masses plotted against rest-frame 3 GHz radio lumi- 
nosity. The small dynamic range of the stellar masses combined with their relatively large 
uncertainties makes the formal determination of any correlation difficult, but a weak correla- 
tion is visually apparent. A proper statistical survival analysis that accounts for the multiple 



^The modeled mass-to-light ratios include the mass of all material that has gone through star-formation. 
It thus represents the integrated star formation history and includes stars, as well as sub-stellar objects 
and compact post-stellar relics such as white dwarfs, neutron stars, and black holes, but not the gas that is 
recycled to the ISM. Stars dominate the rest-frame near-IR SED. 
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upper-limits finds that the Cox Hazard Probabihty is only 0.5 meaning there is only a 50% 
chance that stellar mass correlates with radio luminosity for our sample. The detection of a 
correlation, showing that more powerful radio sources are hosted by more massive galaxies, 
would not be surprising, and is in fact expected if the M — a relation holds to high redshift 
and all HzRGS have similar efficiency rates or Eddington ratios. Previous work, inferring 
stellar masses from the observed near-IR flu x rather than the rest-frame near-IR luminosity, 



has suggested this correlation in the past (lEales et al.l 119971 : iBest et al.l Il998l : IJarvis et al. 



200ll : toe Breuck et all 12001 IWillott et al.ll2003h . 



6.2. Comparisons with More Detailed Modeling 



As mentioned earlier, several HzRGs from our sample have already been the subject 
of detailed studies that included Spitzer observations and determined host stellar masses. 
These studies have generally used more extensive photometry for the modeling, including 
optical and near-IR data, and have used more complicated models to determine the host 
stellar masses. We review the results of more detailed modeling here; in brief, allowing for 
different choices of IMF, the simplistic modeling presented here is consistent with the more 
detailed modeling 



MRC 2104-242 (z = 2.491) was analyzed in detail by IVillar-Martm et al.l (120061 ) who 



found a stellar mass of 5 ± 2 x IO^^Mq, E{B - V) = 0.4, and an age of 1.8 Gyr. This 
result is consistent with the work here where we find a mass of 1.6l^:^ X lO^M©. Our work 
assumes no extinction and Zform = 10 corresponds to an age of 2.1 G y r for the redshift of 
thi s HzRG. Much of the diff erence can be attributed to the ISalpeterl (119551 ) IMF adopted 
by IVillar-Martm et al.l (120061 ) wh ich yield s a rn ass-to-light ratio systematically higher by a 
factor of ~ 40% compared to the iKroupal (120011 ) IMF adopted here. Taking this choice into 
account, the derived masses are within la of each ot her, which is rather encour aging given 
the increased sophistication of the stellar modeling in IVillar-Martm et al.l (120061 ). as well as 
their inclusion of near-IR data from the ground and HST. 



LBDS 53W091 {z = 1.552) was studied by lStern et al.l (120061 ) and, as d iscussed e arlier , 
is (^pr oto-) typical o f the o ld, evolved class of extremely red objects. Using the ISalpeterl (119551 ) 
IMF, Istern et all kood \ find a stellar mass of ~ 3 x lO^^M©, and an age of 3.5 - 4 Gyr. 
Accounting for the difference in IMF, this result is consistent with the masses presented here, 
~ 1.6 X IO^^Mq. The ages of the LBDS galaxies imply a redshift of formation ^form ~ 9, 
similar to our -Zform = 10. 



The full radio to X-ray SED of 4C 23.56 {z = 2.483) is presented in 



De Breuck et al, 
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( 120071 ). the first time such data has been presented for a radio galaxy above z = 2. A full 
SED decomposition of the AGN and stellar emission (such as scatter quasar light, emission 
lines and nebular continuum) is performed. With a Salpeter IMF and no extinction, they 
derive a stellar mass of S.llg i x 10^^ Mq. Again, allowing for the different choice of IMF, 

X 1O"M0. 



this result is consistent within la of the results here, 3.9^n7 x lO^^Mf 



These comparisons give us confidence that though the SED modeling we adopt is some- 
what simplistic, the derived stellar masses appear to be accurate within the uncertainties 
inherent to the enterprise. 



6.3. Systematic Effects on the Stellar Masses 



Deriving stellar masses of high-redshift galaxies {z ^ 1) has become wides pread in recent 



years 



with the advent o f deep near- and mid - IR extra-galactic surve ys {e.g., IShapley et al. 



20031 : iLabbe et al.l l2005l : iRettura et al.l l2006l : iMoustakas et al.l 120071 ). Photometric masses. 



as derived here, are subject to several syste matic uncertainties. H owever, "relative" masses 
can be derived with less uncertainty [e.g., [Dickinson et al.ll2003l ). We now briefly discuss 



some of the systematic effects inherent to the stellar masses derived here. This analysis will 
be particularly important for future comparison of our results with other samples. 



6.3.1. IMF 



The IMF remains one of the largest uncert ainties when m odeling the SED of a galaxy. 



Ther e have been seve ral updates to the classical ISalpeterl (Il955l ) IMF recently {e.g., iKroupa 



200ll : IChabrierll2003l . etc.). Even these, however, are largely based on local stars within our 
own Galaxy and hence suffer uncertainties at the low mass end and for extreme metalicities 
(high and low). We have chosen to use the iKroupal (I20011) IMF here. However, our results 
would not be very d i ffereri t (< 2%) if we used the IChabrieii (120031 ) IMF. Had we used 
the classical ISalpeten (119551 ) IMF, our stellar masses would increase systematically by 40%, 
independent of our choice of Ziorm- 



6.3.2. Formation Redshift 



We have assu med a uniform formation redshift, -Zform = 10, in the current analysis. 

14 for i'-band drop-outs at 2; ~ 6 in the 



Eyles et al.l (120061 ) find a range of Z{o 



7 



GOODS fields for galaxies. These distant galaxies are our best-studied examples of galaxies 
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at 6 ;^ 2; ^ 7 and hence are our best (only) candidates for the parent populations to HzRGs. 
This range of formation redshifts corresponds to 0.16 Gyr earlier or 0.3 Gyr later than z = 10 
when the universe was only ~ 1 Gyr ol d. Assurning th e full range of ^form f rom lEvles et al. 
J2006h changes the mass by ±17% for a lKroupal J200lh IMF and ±8% for a lSalpeteil JlQSsh 
IMF. The change is such that the stellar mass decreases with a later formation redshift and 



vice- versa. 



6.3.3. Extinction 



No internal galactic extinction is considered here. However, ~ 50% of the HzRGs have 
sub- mm detections, implying substantial dust content, though the dust distribution need not 
heavily obscure the old stellar population observed at rest-frame near-IR wavelengths. To 
obtain a 20% under-estimation of the stellar mass frorn the rest-frame H— band luminosity 
requires an Ay > 1.17 assuming a lCalzetti et al.l (119941 ) extinction curve. 



6.3.4. Star-burst Component 

These stellar mass estimates do not include the possible contribution of young stellar 
population which would lead to an over estimation of the stellar mass. If we include a 10% 
by mass starburst component, starting 200 Myr before the observed redshift, we find that 
while the SED at rest-frame near-IR does not change significantly, the mass-to-light ratio 
does decrease by 10 — 20% at i/— band (see Fig. [2]). Redoing our analysis again with this 
younger component results in nearly identical H— band luminosities as before, but a < 20% 
decrease in stellar mass. We note that despite several sources in our sample having sub-mm 
detections, it is not thought that on-going star formation constitutes as much as 10% by 
mass for such systems. 

Deep rest-fram e UV spectropolarimetry of nine optically bright star-forming HzRGs by 



Vernet et al.l (120011 ) finds only one object with a substantial contribution from a population 
of young stars. Hence we conclude that there may be a mild systematic over-estimation 
of the stellar mass, but this uncertainty is likely less significant than other uncertainties 
discussed above, such as the choice of the IMF in the modeling. 



Maraston 



6.3.5. TP- AGB Stars 
(120051 ) has argued that the thermally pulsing asymptotic giant branch (TP- 
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AGB) phase in the evolution of a simple stellar population is not properly conside red in 
other spectral evolution synthesis codes {e.g., PEGASE and iBruzual fc CharlotI l2003l ) . and 
hence analysis using these models over-predicts the near-IR mass-to-light ratio for galaxy 
ages < 1 Gyr. For our choice of ^^form, this age corresponds to galaxies at 3 < z < 4. We 
emphasize that t he PEGASE mode ls do incorporate the TP-AGB phase, just in a different 
manner than the iMarastoru (120051 ) models. If this TP-AGB phase is indeed an issue, it is 
most likely to only affect the highest redshift (and hence youngest) sources in our sample by 
decreasing their stellar mass by < 30%. As we have seen in §6.3.4 an un-accounted sta r-burst 
leads to an overestimation of the stellar mass. Use of the model from iMarastonI (120051 ) would 
hence increase this overestimate to about 30% which is on par with the uncertainty due to 
all the other systematics. 



Discussion and Conclusions 



We have performed a comprehensive infrared survey of 69 HzRGs above z = 1, the 
first time such galaxies have been detected in the rest-frame near-IR and mid-IR. This work 
presents the largest sample of confirmed type 2 quasars to be systematically studied with 
Spitzer. As expected, they predor ninantly res i de in the I RAC color-color sp ace dominated by 
ty pe 1 AGN found ern pirically by lLacy et al.l (120041 ) and lStern et al.l (120051 ) and theoretically 
by ISajina et al.l (120051 ). 



We model the rest-frame near to mid-IR SEDs derived from our broad band IRAC, 
IRS and MIPS imaging. We use a toy model comprised of thermal BBs to represent the 
emission at longer wavelengths presumably due to dust, and the PEGASE code to model 
the stellar component. Our analysis yields total near-IR (if— band) and mid-IR (5 /zm) rest- 
frame luminosities and we deconvolve the stellar and AGN contributions at H— band. The 
scatter of stellar H— band luminosities is lower than that for the total H— band luminosities 
and we find the HzRGs with lower stellar fractions tend to have higher rest-frame H— band 
total luminosities, although this trend is weak. When //—band luminosities are corrected 
for AGN contribution they show no correlation with stellar fraction. Also the stellar fraction 
is not found to correlate with redshift, radio luminosity or mid-IR luminosity. These results 
are consistent with orientation-driven models of AGN unification. 

We find that the H— band stellar luminosities are consistent with stellar masses of 
conip a rable with the mas ses of local massive radio galaxies (iMcCarthyl Il993 



Jarvis et al.l l200ll : IWillott et al.l |2003| ) . This conversion of luminosity to mass assumes a 
dust-free, passively evolving elliptical formed at 2; = 10 with no recent episodes of star- 
formation. We discuss these effects and other systematics and conclude that our conclusions 
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are robust at the level of ~ 10 — 20%. We also note that effects can act to either increase the 
inferred mass (the presence of dust) or to decrease it (recent star-burst activity). So multiple 
effects may cancel out, thus, most likely, simply increasing the scatter in the estimated stellar 
masses. 

The derived stellar masses average lO^^'^M© across z = 1—4 with a scatter less than 50%. 
This scatter is remarkably small, suggesting that the tightness of the observed K — z relation 
at 2 > 1 is due to the intrinsic homogeneity of the HzRG population. The increase in the 
scatter from 2; < 1 to 2; > 1 can largely be att ributed to the -ftT-band being affected by emission 
lines, scattered light from the AGN {e.g., iDe Breuck et al.l 120071 ) and, potentially, young 



stellar populations. The scatter in the observed H— band luminosities can be explained by 
varying contributions from AGN heated dust entering the IRAC bands. 

These masses, while very high, are not uncommon at these redshifts. Much work 
has been done recently to detect massive, M > lO^^Mf^, galaxie s at high redshifts {e.g., 
Labbe et al.l l2005l : iPapovich et al.l l2006l : Ivan Dokkum et al.l |2006| ) . Densities of approxi- 
mately on e M > lO^^MfT) galaxy per square arc-minut e are found in the redshift range 
2 ^ -2 ;^ 3 (IPapovich et al.ll2006l : Ivan Dokkum et al.ll2006l ). Our survey targeted a sample of 
known sources from the literature, derived from multiple radio surveys. As such, we are ill- 
posed to address space densities from the current work. However, we note that our adopted 
definition of a HzRG {z > 1 and Lsghz > lO^^WHz"^) implies that lar ge area survey s such 
as the Faint Images of the Radio Sky at Twenty-centimeters (FIRST; iBeckerl IWhitd ) have 
already detected all HzRGs at 1 < z < 3, though the vast majority lack spectroscopic ob- 
servations thus far. We can therefore constrain the surface density of HzRGs to be less than 
100 per square degree, which is significantly lower than the surface density of optically and 
near-IR selected massive galaxies. This low surface density is most likely due to the short 
lifetimes (10 — 100 Myrs) and low duty cycles of the most powerful radio sources, even if the 
every massive galaxy were to go through a radio-luminous phase. 

We find a slight trend that the galaxies with the higher stellar masses host the more 



powerful radio sources, a trend suggested by other authors (iMcCarthyl Il993l : IJarvis et al. 



2OOII : IWillott et al.ll2003l ). However, a full survival analysis taking proper account of upper 
limits finds this correlation has only a 50% of being correct. The correlation is expected if 
the M — a relation holds to high redshift and HzRGs all accrete at comparable Eddington 
ratios. 

Most of the HzRG mid-IR luminosities are above 10^^ Lq, implying even higher total IR 
luminosities. These LIRG and ULIRG-like luminosities are due to reprocessed radiation from 
gas and dust obscuring a highly accreting AGN. Rest-frame mid-IR luminosity (Aq ~ 5/im) 
may be a good, independent indicator of AGN activity. Comparison with traditional AGN 



- 26 - 



indicators at other wavelengths is currently underway (X-ray — D. Alexander et al., in prep.; 
radio — N. Seymour et al., in prep.). 

We would like to thank the referee for the constructive comments provided on the clar- 
ification of this paper. NS and DS thank ESO for generous hospitality on several occasions 
over the last few years. We thank Dave Prayer and Ranga-Ram Chary for advice on reduc- 
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on 4C 41.17 and 6C 0140+326. This work is based on observations made with the Spitzer 
Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute 
of Technology under a contract with NASA. Support for this work was provided by NASA 
through an award issued by JPL/Caltech. The work by WvB was performed under the 
auspices of the U.S. Department of Energy, National Nuclear Security Administration by 
the University of California, Lawrence Livermore National Laboratory under contract No. 
W-7405-Eng-48. WvB also acknowledges support for radio galaxy studies at UC Merced, 
including the work reported here, with the Spitzer Space Telescope via NASA grants SST 
1264353, 1265551, 1279182 and 1281587. 



A. Notes on Individual Sources 

Here we provide individual comments on selected HzRGs. The SED fitting is presented 
in Pig. [5] (20 sources with MIPS detections) and Pig. [6] (50 without MIPS detection), arranged 
in order of increasing redshift. 



3C65 {z = 1.176): This HzRG was observed in only two bands due to an astrometric 
error in NED which we have since had corrected. 

3C294 {z = 1.786): While this radio galaxy is detected in all IRAC bands, the pho- 
tometry is severely affected by a bright star 8"a to the east by nort h-east. This bright 



star has allowed detailed ada ptive optics imaging of this object (e.g. iQuirrenbach et al. 
200ll : Isteinbring et all [200^ 1 



B3 J2330+3927 {z = 3.086): This HzRG is one of two sources consistent with an AGN 
power law dominating the entire optical to mid-IR emission. We therefore can only 
provide an upper limit to the stellar mass. Although optical a nd near-IR spectroscopy 
( De Breuck et al.ll2003l ) show a classical type 2 AGN spectrum. iPerez- Torres fc De Breuck 
f l2005h :bund that the radio morphology is more consistent with a type 1 AGN. The 
latter would be in agreement with the high contribution from hot AGN light seen in 
the near- to mid-IR SED. 
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• MRC 0316—257 {z = 3.131): Due to an error at the time of proposal sub mission, the 
source observed with IRAC was MRC 0316— 257B ( Le Fevre et ahl 1996 ). a spectro- 
scopically confirmed z = 3.118 AGN, 71" from MRC 0316—257. We used the deeper 
Cycle 1 GO data (Program ID 3482, PI van Breugel), but no IRS 16 fim data is avail- 
able for the z > 2 radio galaxy. 

• TX J1908+7220 {z = 3.53): This HzRG is the second of two sources consistent with 
an AGN power law dominating the entire optical to mid-IR emission. We therefore 
can only provide an upper limit to the stellar mass. 



TN J1338-1942 



4. 11): T his HzRG is the fourth highest redshift source in our 



sample (IDe Breuck et al.lll999l ). Although the 3.6 /xm IRAC band is likely to be con- 
taminated by Ha, the remaining bands are consistent with a stellar SED. The low 
8.0 /im flux ensures that any hot dust contamination will be small. 

8C 1435-1-635 (z = 4.25): This HzRG is the third highest redshift source in our sample 



( iLacy et al.lll994l ). The stellar mass estimate is based on the 4.5 /xm IRAC detection 
only, as the 3.6 /im band is likely to be contaminated by Ha emission. 

6C 0140+326 {z = 4.41): This HzRG is the second most distant galaxy and is masked 
by a bright, foreground galaxy in IRAC channels one and two, thus rendering source 
extraction at these wavelengths infeasible. 

TN J0924— 2201 (z = 5.1 9): This HzRG is the most distant radio galaxy in our sample 



( Ivan Breugel et al.lll999l ). It is detected only at 3.6 and 4.5 /im. Deeper IRAC data 
has been obtained (Program ID 20749, PI Zheng), but these data were not public yet 
at the time of writing (December 2006). 
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Table 1. Spitzer HzRG Sample and Exposure Times per Instrument 





RA 


Dec. 




log(I/3GHz) 


/f-band 


IRAC 


IRS 


MIPS 


HzRG 


(,T2000) 


(.T2000) 


R.cdshift 


(W Hz~ ) 


(Vega) 


(s) 


(s) 


(s) 


6C 0032+412 


00:34:53.1 


41:31:31.50 


3.670 


27.75 


19.10 


120 


122 


134, 420, 881 


MRC 0037-258 


00:39:56.4 


—25:34:31.01 


1.100 


27.07 


17.10 


120 








U.L.UJ-..LO.f 


AQ' 9 9Q 


1.173 


26.68 


17.60 


120 




134 420 881 


MRC 0114-211 


01:16:51.4 


—20:52:06.71 


1.410 


28.39 


18.50 


120 






TN J0121+1320 


01:21:42.7 


13:20:58.00 


3.516 


27.41 


18.70 


120 


122 




6C*0132+330 


01:35:30.4 


33:17:00.82 


1.710 


26.60 


18.80 


120 






6C 0140+326 


01:43:43.8 


32:53:49.31 


4.413 


27.89 


19.90 


5000 


122 


267, 671, 2643 


MRC 0152-209 


01:54:55.8 


— 20:40:26.30 


1.920 


27.77 


17.90 


120 






MRC 0156-252 


01:58:33.6 


-24:59:31.10 


2.016 


27.79 


16.10 


120 


122 




TN J0205+2242 


02:05:10.7 


22:42:50.40 


3.506 


27.43 


18.70 


120 


122 




MRC 0211-256 


02:13:30.5 


— 25:25:21.00 


1.300 


26.22 




120 


122 




3C65 


02:23:43.2 


40:00:52.40 


1.176 


28.06 


16.80 


120 






MRC 0251-273 


02:53:16.7 


— 27:09:13.03 


3.160 


28.09 


18.50 


120 


122 




MRC 0316-257 


03:18:12.0 


— 25:35:11.0(3 


3. 13(3 


28.26 


18.10 


120 


122 




MRC 0324-228 


03:27:04.4 


— 22:39:42.60 


1.894 


27.81 


18.80 


120 






MRC 0350-279 


03:52:51.6 


—27:49:22.61 


1.900 


27.61 


19.00 


120 




134 420 881 


MRC 0406-244 


04:08:51.5 


— 24:18:16.39 


2.427 


28.11 


17.40 


120 


122 


134, 420, 881 


4C60.07 


05:12:54.8 


60:30:52.01 


3.788 


27.91 


19.20 


120 


122 


134, 420, 881 




05:30:25.2 


54:54: 22. 0(3 


2. 575 


28.27 




120 


122 


134 420 881 


WN J0617+5012 


(36:17:39.4 


5(3: 12: 55.4(3 


3.153 


26.97 


19.70 


120 


122 




4C41.17 


06:50:52.1 


41:3(3:31.(3(3 


3.792 


28.17 


19.10 


5000 


122 


267 671 2643 


WN J0747+3654 


07:47:29.4 


36:54:38.09 


2.992 


27.02 


20.00 


120 


122 




6CE0820+3642 


08:23:48.1 


36:32:46.42 


1.860 


27.49 


18.20 


120 






IT*^^; nS9R-l-1<3'^ 

*J OO UO^OT^ J- y o 


L/O.OL/.UO.'i 


1 Q-1 "^-i fi nn 


2.572 


27.47 


18.20 


120 


122 




5C7.269 


08:28:38.8 


25:28:27.10 


2.218 


27.08 


18.90 


120 


122 




6CE0901+3551 


09:04:32.4 


35:39:03.23 


1.910 


27.47 


18.10 


120 






6CE0905+3955 


09:08:16.9 


39:43:26.(30 


1 .883 


27.49 


18.30 


120 






B2 0902+34 


09:05:30.1 


34:(37:56.89 


3.395 


28.27 


19.90 


1200 


122 


2557, 2696, 3556 


TN J0924-2201 


09:24:19.9 


— 22:01:41.00 


5.195 


27.83 


19.70 


120 


122 


134 420 881 


fin flQ^fl-l-^RQ 




'\fi.■A^ - 1 4 

00.*±X.iJL/.-L*± 


2.395 


27.79 


19.50 


120 


122 








— 24:28:49.65 


2.923 


27.95 


19.20 


120 


122 


134, 420, 881 


3C239 


10:11:45.4 


46:28:19.75 


1.781 


28.19 


17.80 


120 






MRC 1017-220 


10:19:49.0 


— 22:19:58.03 


1.768 


28.11 


17.40 


120 






MG 1019 + 0534 


10:19:33.4 


5:34:34.80 


2.765 


28.13 


19.10 


120 


122 




WN Jllln + 5016 


11:15:06.9 


50:16:23.92 


2.540 


26.87 


19.20 


120 


122 




3C257 


11:23:09.2 


5:30:19.47 


2.474 


28.62 


17.80 


120 


122 




WN J1123+3141 


11:23:55.9 


31:41:26.14 


3.217 


27.42 


17.40 


120 


122 




PKS 1138-262 


11:40:48.6 


— 26:29:08.50 


2.156 


28.14 


16.10 


3000 


122 


9000, -, - 


3C266 


11:45:43.4 


49:46:08.24 


1.275 


27.80 


17.10 


120 






6C 1232+39 


12:35:04.8 


39:25:38.91 


3.220 


28.01 


18.00 


120 


122 




USS 1243 + t)36 


12:45:38.4 


3:23:20.70 


3.570 


28.25 


19.20 


120 


122 




TN J1338-1942 


13:38:26.0 


— 19:42:31.(3(3 


4. 1 1(3 


27.90 


19.60 


5000 


122 




4C24.28 


13:48:14.8 


24:15:52.00 


2.879 


28.25 




120 


122 




3C294.0 


14:06:44.0 


34:11:25.00 


1.786 


28.12 


17.90 


120 






USS 1410-001 


14:13:15.1 


0:22:59.70 


2.363 


27.69 




120 


122 




SC 1435+635 


14:36:37.1 


63:19:14.00 


4.250 


28.55 


19.40 


120 


122 


134 420 881 




16:01 : 17.3 




2.527 


28.00 




120 


122 




USS 1707+105 


17:10:06.5 


10:31:06.00 


2.349 


27.78 




120 


122 




LBDS 53w002 


17: 14:14.7 


50:15:29.70 


2.393 


27.02 


18.90 


3300 


122 


134, 420, 881 


LBDS 53w091 


17:22:32.7 


50:06:01.94 


1.552 


26.29 


18.70 


900 




633, 1311, 2643 


3C366.0 


17:24:19.0 


50:57:40.30 


1.079 


27.65 


16.80 


120 




134, 420, 881 


7C 1751+6809 


17:50:49.9 


68:08:25.93 


1.540 


27.01 


18.20 


120 




134, 420, 881 


7C 1756+6520 


17:57:05.4 


65:19:53.11 


1.480 


27.00 


18.90 


120 




134, 420, 881 


3C36S.0 


18:05:06.3 


11:01:33.00 


1.132 


27.63 


17.20 


120 




134, 42(3. 881 


7C 1805+6332 


18:05:56.9 


63:33:13.14 


1.840 


27.12 


18.80 


120 




134, 420, 881 


4C40.36 


18:10:55.7 


40:45:24.01 


2.265 


27.94 


17.80 


120 


122 


134, 420, 881 


TX J1908+7220 


19:08:23.7 


72:20:11.82 


3.530 


28.15 


16.50 


120 


122 


134, 420, 881 


WN J1911+6342 


19:11:49.6 


63:42:09.60 


3.590 


27.03 


19.90 


120 


122 


134, 420, 881 


TN J2007-1316 


20:07:53.3 


-13:16:43.62 


3.840 


27.79 


18.80 


120 


122 




MRC 2025-218 


20:27:59.5 


-21:40:56.90 


2.630 


27.96 


18.50 


120 


122 




MRC 2048-272 


20:51:03.6 


-27: (33: (32. 53 


2.(36(3 


27.85 


18.30 


120 


122 




MRC 2104-242 


21:06:58.1 


-24:05:11.00 


2.491 


27.88 




120 


122 




4C23.66 


21:07:14.8 


23:31:45.00 


2.483 


28.26 


19.73 


120 


122 


134, 420, 881 


MG 2144+1928 


21:44:07.5 


19:29:14.60 


3.592 


28.27 


19.10 


120 


122 




USS 2202+128 


22:05:14.1 


13:05:33.50 


2.706 


27.75 


18.40 


120 


122 




MRC 2224-273 


22:27:43.3 


-27:05:01.71 


1.679 


27.39 


18.50 


120 
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Table 1 — Continued 





RA 


Dec. 




log(i3GHz) 


/f-band" 


IRAC 


IRS 


MIPS 


HzRG 


(J2000) 


(J2000) 


Redshift 


(W 


(Vega) 


(s) 


(s) 




B3 ,72330 + 3927 


23:30:24.9 


39:27:12.02 


3.086 


27.59 


18.80 


120 


122 


134, 420, 881 


4C28.58 


23:51:59.2 


29:10:28.99 


2.891 


27.89 


18.70 


120 


122 




3C470 


23:58:35.3 


44:04:38.87 


1.653 


28.24 


18.50 


120 




134, 420, 881 



^/T-band magnitude within projected 64kpc radius 
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Table 2. Spitzer Photometric Parameters 





Wavelength 


Pixel Scale 


Conv. 


Aperture" 


Annulus** 


Aperture 


Filter 


(urn) 


(7pix) 


to /iJy 


(") 


(") 


Correction 


IRACl 


3.6 


0.61 


23.5 


3.5 


9.76 


1.48 










7.0 


9.76 


1.13 


IRAC2 


4.5 


0.61 


23.5 


3.5 


9.76 


1.60 










7.0 


9.76 


1.168 


IRAC3 


5.8 


0.61 


23.5 


3.5 


9.76 


2.76 










7.0 


9.76 


1.155 


IRAC4 


8.0 


0.61 


23.5 


3.5 


9.76 


2.13 










7.0 


9.76 


1.316 


IRS-blue 


16 


1.2 


0.59r 


6 


6-10 




MIPSl 


24 


1.25 


36.71 


13 


20-32 


1.167 


MIPS2 


70 


4.95 


576 


35 


39-65 


1.211*^ 


MIPS3 


160 


8 


1500 


50 


75-125 


1.445 



^Diameter in which the photometry is measured. 

'^Inner and outer diameter of annulus used to determine the background sky level, 
except for the IRAC bands where we provide the size of the background mesh used by 
SExtractor. 

"^Including aperture correction. 

'^For a power-law SED, oc 
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Table 3. Spitzer Photometry for HzRGs (/xJy) 



HzRG 




/4.5/xm 




/s.O/^m 




/24^m 




/l60^m 


6C 0032 + 412 


15.4 ± 2.1 


33.8 ± 3.7 


57.6 ± 8.9 


98.3 ± 10.0 


205.0 ± 34.9 


457.0 ± 59.7 


< 3950 


< 96000 


MR C0037-258 


221.0 ± 22.0 


248.0 ± 25.0 


286.0 ± 29.0 


518.0 ± 52.0 


— 


— 


— 


— 


60*0058+496 


82.0 ± 8.3 


86.7 ± 8.8 


93.2 ± 9.5 


309.0 ± 31.0 


— 


1460.0 ± 87.3 


18900 ± 1962 


< 94600 


MRC 0114-211 


87.3 ± 8.9 


117.0 ± 12.0 


157.0 ± 16.0 


398.0 ± 40.0 


— 


— 


— 


— 


TN J0121+1320 


10.5 ± 1.8 


14.4 ± 2.1 


< 32.5 


< 40.3 


< 136.0 


— 


— 


— 


60*0132+330 


31.6 ± 3.6 


41.3 ± 4.5 


51.9 ± 12.4 


109.0 ± 11.0 


— 


— 


— 


— 


6C 0140+326 


< 623.0 


< 450.0 


34.9 ± 2.7 


36.3 ± 3.4 


194.0 ± 40.0 


2190.0 ± 110.5 


< 2580 


< 20700 


MRC 0152-209 


108.0 ± 11.0 


165.0 ± 17.0 


215.0 ± 22.0 


415.0 ± 42.0 


— 


— 


— 


— 


MRC 0156-252 


291.0 ± 29.0 


405.0 ± 41.0 


717.0 ± 72.0 


1125.0 ± 113.0 


1430.0 ± 116.6 


— 


— 


— 


TN J0206+2242 


7.6 ± 1.7 


4.7 ± 2.3 


< 29.3 


< 37.4 


< 122.0 


— 


— 


— 


MRC 0211-266 


166.0 ± 17.0 


197.0 ± 20.0 


222.0 ± 22.0 


278.0 ± 28.0 


330.0 ± 42.6 


- 


- 


- 


3C65 


161.0 ± 32.0 


- 


320.0 ± 66.0 


- 


- 


- 


- 


- 


MRC 0251-273 


11.5 ± 1.6 


12.7 ± 1.8 


< 29.1 


< 28.8 


97.2 ± 34.2 


- 


- 


- 


MRC 0316-257 


19.3 ± 2.1 


20.1 ± 2.1 


19.5 ± 2.4 


38.1 ± 4.1 


< 91.5 


— 


— 


— 


MRC 0324-228 


39.4 ± 4.2 


39.7 ± 4.3 


61.1 ± 8.6 


89.9 ± 9.9 


— 


— 


— 


— 


MRC 0350-279 


23.6 ± 2.6 


40.6 ± 4.2 


82.2 ± 27.6 


79.3 ± 26.5 


- 


321.0 ± 60.7 


< 1040 


< 88800 


MRC 0406-244 


40.4 ± 4.3 


43.3 ± 4.6 


< 51.6 


63.5 ± 14.5 


242.0 ± 34.6 


1420.0 ± 83.7 


24700 ± 2306 


< 47700 


4C60.07 


27.3 ± 3.1 


33.6 ± 3.7 


35.1 ± 8.4 


37.2 ± 9.5 


175.0 ± 33.0 


1340.0 ± 83.7 


< 3750 


< 64600 


PKS 0529-549 


46.6 ± 4.9 


52.9 ± 5.5 


62.7 ± 8.6 


72.2 ± 9.1 


134.0 ± 33.7 


942.0 ± 71.0 


< 4110 


< 74100 


WN J0617+6012 


3.6 ± 1.0 


6.5 ± 1.2 


< 48.3 


< 54.5 


< 93.4 


- 


- 


- 


4C41.17 


23.4 ± 2.4 


27.5 ± 2.8 


35.6 ± 3.7 


36.5 ± 3.5 


< 99.3 


370.0 ± 43.3 


< 3210 


< 26300 


WN J0747+3654 


19.1 ± 2.4 


25.3 ± 3.0 


< 29.9 


44.9 ± 11.5 


< 111.0 


- 


- 


- 


6CE0820+3642 


79.2 ± 8.1 


81.9 ± 8.4 


82.0 ± 8.4 


68.0 ± 7.0 


— 


— 


— 


— 


USS 0828+193 


61.7 ± 6.9 


133.0 ± 13.0 


201.0 ± 21.0 


687.0 ± 74.0 


1270.0 ± 112.8 


— 


— 


— 


5C7.269 


41.0 ± 4.5 


49.5 ± 5.3 


57.8 ± 6.1 


< 40.1 


< 186.0 


- 


- 


- 


6CE0901+3551 


37.2 ± 4.1 


46.5 ± 5.0 


52.8 ± 10.7 


69.8 ± 12.4 


- 


- 


- 


- 


6CE0905+3955 


51.8 ± 5.4 


60.1 ± 6.2 


96.8 ± 9.8 


146.0 ± 14.0 


— 


— 


- 


— 


B2 0902+34 


6.4 ± 0.8 


9.9 ± 1.1 


11.0 ± 2.5 


41.3 ± 2.3 


144.0 ± 39.6 


323.0 ± 18.8 


— 


< 18500 


TN J0924-2201 


11.3 ± 1.8 


11.4 ± 1.9 


< 30.6 


< 33.0 


< 108.0 


< 160.0 


< 3330 


< 52900 


6C 0930+389 


30.7 ± 3.4 


32.2 ± 3.6 


37.3 ± 9.2 


< 32.6 


< 97.2 


— 


— 


— 


USS 0943-242 


21.5 ± 2.6 


28.4 ± 3.2 


< 30.9 


26.8 ± 11.7 


109.0 ± 34.1 


493.0 ± 57.5 


< 3390 


< 50900 


3C239 


96.4 ± 9.8 


111.0 ± 11.0 


130.0 ± 12.0 


142.0 ± 14.0 


— 


— 


— 


— 


MRC 1017-220 


119.0 ± 12.0 


179.0 ± 18.0 


273.0 ± 27.0 


360.0 ± 36.0 




— 


— 


— 


MG 1019+0534 


25.6 ± 2.9 


19.5 ± 3.8 


< 36.4 


< 42.9 


■ ± 52.5 


— 


— 


— 


WN J1115+5016 


7.8 ± 2.1 


9.5 ± 3.0 


< 54.8 


< 61.1 


< 88.0 


- 


- 


- 


3C257 


85.0 ± 8.7 


111.0 ± 11.0 


194.0 ± 19.0 


322.0 ± 33.0 


569.0 ± 74.6 


- 


- 


- 


WN J1123+3141 


48.2 ± 5.0 


74.4 ± 7.6 


92.7 ± 9.4 


182.6 ± 18.4 


742.0 ± 70.4 


- 


- 


- 


PKS 1138-262 


318.0 ± 32.0 


497.0 ± 50.0 


887.0 ± 89.0 


1500.0 ± 150.0 


2000.0 ± 160.3 


3890.0 ± 176.2 


— 


— 


3C266 


67.9 ± 7.0 


73.1 ± 7.5 


45.1 ± 4.7 


102.6 ± 10.4 


- 


- 


- 


- 


6C 1232+39 


33.3 ± 3.6 


41.8 ± 4.4 


52.0 ± 9.2 


75.2 ± 10.8 


208.0 ± 35.9 


- 


- 


- 


USS 1243+036 


22.0 ± 2.5 


21.5 ± 3.0 


< 48.0 


< 62.9 


< 163.0 


— 


— 


— 


TN J1338-1942 


17.8 ± 1.9 


10.7 ± 1.2 


14.3 ± 2.8 


9.9 ± 3.4 


< 119.0 


- 


- 


- 


4C24.28 


16.5 ± 2.2 


27.3 ± 3.1 


43.3 ± 8.3 


102.0 ± 10.0 


446.0 ± 50.8 


- 


- 


- 


3C294.0 


< 93.0 


< 103.0 


68.0 ± 16.8 


66.6 ± 20.6 


- 


- 


- 


- 


USS 1410-001 


50.6 ± 5.3 


79.0 ± 8.1 


166.0 ± 17.0 


240.0 ± 24.0 


481.0 ± 58.6 


— 


— 


— 


8C1436+635 


18.6 ± 2.1 


14.7 ± 2.9 


< 45.5 


< 50.3 


< 86.7 


< 165.0 


< 4130 


< 33900 


USS 1558-003 


78.8 ± 8.1 


101.0 ± 10.3 


105.0 ± 10.5 


233.0 ± 23.4 


398.0 ± 52.2 


— 


- 


- 


USS 1707+105 


22.1 ± 2.7 


30.1 ± 3.4 


22.6 ± 9.1 


< 33.3 


< 123.0 








LBDS 53w002 


32.0 ± 3.3 


44.0 ± 4.5 


49.9 ± 5.2 


103.0 ± 11.0 


177.0 ± 33.9 


648.0 ± 61.4 


< 4300 


< 65100 


LBDS 53w091 


43.3 ± 6.3 


51.4 ± 5.3 


23.9 ± 6.1 


26.6 ± 6.4 




< 45.9 


< 1610 


< 23900 


3C356.0 


108.0 ± 11.0 


110.0 ± 11.0 


122.0 ± 14.0 


434.0 ± 47.0 




4060.0 ± 192.3 


< 4400 


< 70200 


7C 1761+6809 


46.6 ± 4.9 


50.8 ± 5.3 


< 40.9 


36.5 ± 16.0 




351.0 ± 50.5 


< 3600 


< 51600 


7C 1766+6520 


39.6 ± 4.2 


46.9 ± 6.0 


34.7 ± 7.7 


42.6 ± 8.6 




444.0 ± 51.1 


< 5900 


< 151000 


3C368.0 


126.0 ± 13.0 


112.0 ± 11.0 


112.0 ± 11.0 


210.0 ± 21.0 




3250.0 ± 166.7 


28800 ± 2710 


< 39000 


7C 1805+6332 


28.4 ± 3.6 


42.1 ± 4.4 


51.4 ± 5.4 


95.6 ± 17.1 




648.0 ± 56.2 


< 4330 


< 111000 


4C40.36 


36.5 ± 3.9 


41.3 ± 4.3 


45.4 ± 12.9 


26.3 ± 10.3 


115.0 ± 29.4 


520.0 ± 55.2 


< 4750 


< 64500 


TX J1908+7220 


200.0 ± 20.0 


229.0 ± 23.0 


241.0 ± 25.0 


480.0 ± 48.0 


841.0 ± 70.6 


1910.0 ± 98.9 


16200 ± 1905 


< 63300 


WN J1911+6342 


9.6 ± 1.8 


8.7 ± 1.8 


< 26.4 


< 26.0 


< 75.3 


< 142.0 


< 3710 


< 57100 


TN J2007-1316 


56.2 ± 5.9 


54.4 ± 5.8 


41.1 ± 10.7 


121.3 ± 12.3 


244.0 ± 53.5 








MRC 2025-218 


68.4 ± 7.1 


77.1 ± 8.0 


86.8 ± 10.9 


126.8 ± 12.9 


190.0 ± 50.2 








MRC 2048-272 


59.5 ± 6.2 


72.6 ± 7.5 


78.3 ± 10.4 


38.4 ± 13.1 


< 141.0 








MRC 2104-242 


28.1 ± 3.3 


29.7 ± 3.5 


32.8 ± 10.0 


< 36.3 


121.0 ± 47.0 








4C23.66 


61.1 ± 6.4 


86.2 ± 8.8 


126.9 ± 12.8 


423.7 ± 42.5 


1610.0 ± 129.8 


4390.0 ± 203.8 


30300 ± 2958 


< 70600 


MG 2144+1928 


22.1 ± 2.7 


18.3 ± 2.4 


< 25.5 


< 30.3 


< 93.9 








USS 2202+128 


60.4 ± 6.3 


95.8 ± 9.8 


120.1 ± 12.1 


106.6 ± 10.8 


307.0 ± 43.2 








MRC 2224-273 


61.6 ± 6.4 


86.1 ± 8.8 


98.4 ± 10.0 


203.9 ± 20.5 










B3 J2330+3927 


99.6 ± 10.1 


143.0 ± 14.0 


160.0 ± 16.0 


474.0 ± 47.0 


1040.0 ± 88.0 


2210.0 ± 113.2 


< 4670 


< 64300 
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Tablc 3 — Continued 



HzRG 


/3.6/im 


/4.5^m 


/5.8/im 


/s.O/im 


/l6;Lim 


/24/Lim 


/70/Lim 


/l60/[xm 


4C28.B8 


31.6 ± 3.6 


36.0 ± 3.9 


< 41.7 


40.9 ± 4.4 


228.0 ± 31.5 








3C470 


49.5 ± 10.4 


75.2 ± 11.8 


70.9 ± 10.4 


266.0 ± 30.0 




2660.0 ± 133.5 


< 5570 


< 102000 
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Table 4. Results of SED fitting for HzRGs with MIPS detections 





'\nrc( r tot It \ 


T Stel It \ 


J stel 


log^ikZ /ikzQ j 


T 1 T \ 


6C 0032+412 


12.13 


11.47+°-}^ 

1 _o,12 


0.22 


11.28+nlo 


11.97+no5 
J. j-.^ 1 _o.24 


6C*0058+495 


11.16 


^—0.10 


0.99 


11.26+°-l^ 

"—0.10 


11.40+°-}? 

"—0.11 


MRC 0350-279 


11.29 


10.89+°-?? 

— U.i7 


0.40 


10.85+°-?? 

— U.i7 


10.91+°-^° 


MRC 0406-244 


11.50 


11 49+°-!i^ 


0.98 


11 38+^-!? 


11 75+^-!n 


4C60.07 


11.67 


11.65+°-|^ 

— U.IU 


0.95 


11.44+°-!;^ 

— u.lu 


12.50+^-^f 


PKS 0529-549 


11.67 


11 60+^1^ 


0.86 


11 46+^-1?, 

-^-^■^'-'— 0.10 


11 60+^-}^ 

J. j..uu_Q 10 


4C41.17 


11.60 


11.0U_o 08 


0.99 


11-391°:°^ 


12.061°:^^ 


B2 0902+34 


11.52 


10.97t|]:Jt 


0.28 


io.8it°:lt 


11.801°:!^ 


USS 0943-242 


11.40 




0.99 


11 22+°-i5 

±i.ZZ_Q 


1 1 67+0-35 

i±.U ( _Q 20 


PKS 1138-262 


12.59 


< 12.18 


< 0.39 


< 12.11 


1 2 '^Q+0-35 


LBDS 53w002 


11.55 


11.381°:- 


0.68 




11.551°:- 


3C356.0 


11.28 


11.281°:!^ 


1.00 


1±.0»_Q IQ 


11.391°:°^, 


7C 1751+6809 


11.14 


11-121^:}^ 


0.95 


11 14+0-13 


10.391°:?° 


7C 1756+6520 


11.12 


11 12+°-i2 


1.00 


11 1^^+0-12 
ii.iO_o.io 


10.80+g-?? 


3C368.0 


11.32 


11 32+°-ii 


0.99 


ii.43i°:J^ 


11.281°:^° 


7C 1805+6332 


11.22 


11.111°:}^ 


0.78 


11-071°:}? 


11 19+0-17 


4C40.36 


11.38 


ii.38i°:}3 


1.00 


11 29+°-i3 
ii.zy_Q iQ 


11 21+°-23 
ii.Zi_o.l5 


TX J1908+7220 


12.63 


< 12.44 


< 0.64 


< 12.27 


12.60l°:i? 


4C23.56 


11.91 


n 71 +0-12 


0.63 


11 ^^9+0-12 


19 4Q+0-10 

iz.^y_oo9 


B3 J2330+3927 


12.30 


< 12.11 


< 0.64 


< 11.94 


12.581°:?^ 


3C470 


11.32 


11.301°:?^ 


0.95 


11.301°:?^ 


11.801°:^^ 



-41 - 



Table 5. Results of SED fitting for HzRGs with no MIPS detections 



HzRG logflff /I : ) log(lCi ) log(Lr,fnn/L . ) 



MRC 0037-258 


< 11.56 


< 11.70 


- 


MRC 0114-211 


< 11.37 


< 11.40 


- 


TN J0121+1320 


11.21 


11.00 


- 


60*0132+330 


11.11 


11.10 


- 


6C 0140+326 


< 11.78 


< 11.40 


- 


MRC 0152-209 


< 11.72 


< 11.70 


- 


MRC 0156-252 


< 12.23 


< 12.20 


12.21 


TN J0205+2242 


< 10.96 


< 10.80 


- 


MRC 0211-256 


< 11.54 


< 11.60 


11.23 


3C65 


11.44 


11.50 


- 


MRC 0251-273 


11.14 


11.00 


- 


MRC 0316-257 


11.39 


11.20 


- 


MRC 0324-228 


11.21 


11.20 


- 


WN J0617+5012 


10.74 


10.60 


- 


WN J0747+3654 


11.34 


11.20 


- 


6CE0820+3642 


11.54 


11.50 


- 


USS 0828+193 


< 11.72 


< 11.60 


12.34 


5C7.269 


11.48 


11.40 


- 


6CE0901+3551 


11.32 


11.30 


- 


6CE0905+3955 


11.40 


11.40 


- 


TN J0924-2201 


11.48 


11.10 


< 11.79 


6C 0930+389 


11.40 


11.30 


- 


3C239 


11.60 


11.60 


- 


MRC 1017-220 


< 11.74 


< 11.70 


- 


MG 1019+0534 


< 11.34 


< 11.20 


11.50 


WN J1115+5016 


< 10.86 


< 10.70 


- 


3C257 


< 11.83 


< 11.70 


11.96 


WN J1123+3141 


< 11.81 


< 11.60 




3C266 


11.22 


11.30 




6C 1232+39 


< 11.66 


< 11.50 




USS 1243+036 


< 11.47 


< 11.30 




TN J1338-1942 


< 11.26 


< 11.00 




4C24.28 


< 11.20 


< 11.00 




3C294.0 


< 11.50 


< 11.50 




USS 1410-001 


< 11.53 


< 11.40 


11.86 


8C 1435+635 


< 11.45 


< 11.10 


< 11.63 


USS 1558-003 


11.85 


11.70 


11.83 


USS 1707+105 


11.28 


11.20 
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Table 5 — Continued 



HzRG 


log(L*H°Vi^0) 


log(M^*|^/Mo) 


log(L5^m/-^0) 


LBDS 53w091 


11.20 


11.20 


< 10.88 


WN J1911+6342 


< 11.13 


< 10.90 


< 11.33 


TN J2007-1316 


< 11.79 


< 11.60 




MRC 2025-218 


11.80 


11.60 


11.53 


MRC 2048-272 


< 11.55 


< 11.50 




MRC 2104-242 


11.34 


11.20 


11.30 


MG 2144+1928 


< 11.48 


< 11.30 




USS 2202+128 


< 11.85 


< 11.70 


11.77 


MRC 2224-273 


< 11.34 


< 11.30 




4C28.58 


< 11.56 


< 11.40 


11.68 
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Fig. 1. — Hubble K — z diagram of radio- loud (pluses and filled circles) and radio-quiet 
galaxies (points). The radio-loud sample is from De Breuck et al. (2002). The radio quiet 

samples are from the Hawaii survey (Songaila ct al. 1994) and the Hubble Deep Field - North 
(Dickinson et al. 2003). The filled circles indicate the 69 2; > 1 radio galaxies which comprise 
our Spitzer sample. The three HzRGs falling below the obvious correlation are MRC 0156- 
252 {z = 2.016), PKS 1138-262 {z = 2.156)and TX J1908+7220 {z = 3.530). MRC 0156-252 
is known to have strong Ha in the i^-band (Eales & Rawlings 1996) and PKS 1138-262 has 
also recently been confirmed to have strong Ha in the X-band (Nesvadba et al. 2006). TX 
J1908-I-7220 has a strong unresolved AGN component (see notes in individual sources) and 
seems closer to being an unobscured quasar than a radio galaxy. 
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Fig. 2. — The light-to- mass ratios (in solar units) as a function of age for the adopted 
PEGASE.2 elliptical galaxy model for different rest-frame optical (SDSS m', g\ r') and near- 
IR (2MASS H) passbands. The light-to-mass ratio for the H— band filter, near the peak 
of stellar emission, varies much less as a function of age than the shorter wavelength filters. 
Additionally, the optical passbands are more susceptible to (i) dust absorption — the relative 
effects of one magnitude of visual extinction in each of the passbands is illustrated by the 
horizontal lines; in units of A\, iJ— band is more than four times less sensitive to dust 
extinction than r', (ii) uncertainty in the age of the stellar population, and (iii) the addition 
of light from a second, younger population. For a formation redshift of Zform = 10, the 
corresponding ages are marked near the bottom of the figure for different redshifts. For 
^form = 14 or 7, the tracks shift 0.2Gyrs to the left or O.SGyrs to the right, respectively. 
As the iJ— band track is fairly flat, this uncertain formation redshift has little effect on our 
derived light-to-mass ratio. 
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Fig. 3. — Radio luminosity at rest-frame 3 GHz plotted against redshift for 263 z > 1 radio 
galaxies from the literature (pluses). Filled symbols indicate sources in our iS^iteer program, 
carefully chosen to uniformly cover 1.5 — 2 orders of magnitude in radio power at each redshift 
1 < z < 4. The entire Spitzer sample was observed with IRAC. Radio galaxies at 2; > 2 were 
imaged at 16 /xm with IRS. Radio galaxies in "low" Galactic infrared background (see §3) 
are indicated as sohd squares and were observed with MIPS. 
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-1 1 2 3 -0.4 -0.2 0.2 0.4 0.6 



[5.8] - [B.O] (Vega) log(S5.e/S3 ,) 

Fig. 4. — Mid-IR colors of HzRGs. Solid lines indicate empirical criteria that separate (type 
1) active galaxies from Galactic stars and normal galaxies. The left-hand panel shows the 
criteria from Stern et al. (2005); the right-hand panel shows the criteria originally from Lacy 
et al. (2004), but updated in Lacy et al. (2006). Note that the axes have been selected to 
facilitate comparison with the original works. Crosses indicate typical uncertainties on the 
mid-IR photometry. HzRGs are the quintessential type 2, luminous AGN. As expected, they 
generally reside with the empirical AGN wedges. Tracks from redshift to 2 for a late type 
star-burst galaxy (M82; dashed line) and an early type galaxy (NGC4429; dotted line) are 
included. These template SEDs are from Devriendt et al. (1999) and circles mark the z = 
starting points and asterisks indicate z = 1 and 2; = 2 for the tracks. The tracks indicate 
that neither galaxy would be selected as an AGN candidate if located below z = 2. 
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Fig. 5. — SED fitting of 21 sources with both IRAC observations and MIPS detections. 
The name of each HzRG and its redshift are marked at the top of the panel. Observed 
wavelength is along the bottom axis and rest-frame wavelength is along the top axis. Data 
with error-bars are presented as filled circles: open circles are not used in the fitting for 
reasons described in the text. Downward arrows indicate upper limits. The solid dark line 
indicates the total best-fit SED. The stellar and dust components are indicated by dashed 
and dotted lines respectively. 
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Fig. 6. — Galaxy SED fits to determine stellar masses or upper-limits. The name of each 
HzRG and its redshift are marked at the top of the panel. Observed wavelength is along 
the bottom axis and rest-frame wavelength is along the top axis. Data with error-bars are 
presented as filled circles; open circles not used in the fitting are described in the text. 
Downward arrows indicate upper limits. The dashed line indicates the total best-fit SED. 
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Fig. 7. — Fraction of the HzRG rest-frame H— band luminosity which is modeled as stellar, 
plotted as a function of total H— band luminosity and stellar H— band luminosity. The 
stellar fraction appears independent of the stellar near-IR luminosity but is inversely related 
to the total near-lR luminosity. This implies that the underlying hosts constitute a homoge- 
neous, similar stellar mass population, but with varying AGN contributions to their observed 
total near-IR luminosities. 
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Fig. 8. — Fraction of the HzRG rest-frame H— band luminosity which is modeled as stellar, 
plotted as a function of (i) redshift, (ii) rest-frame 3 GHz radio power, and (iii) rest-frame 
5/im luminosity. Only the 21 HzRGs with MIPS detections are plotted. No correlations are 
evident. 
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Fig. 9. — Rest-frame H— band stellar luminosity versus redshift for the SpitzerYiiRG sample, 
derived from the best-fit models to the multi-band photometry. Solid circles indicate those 
luminosities derived from HzRGs with MIPS detections, whilst open circles are luminosities 
derived from HzRGs without MIPS detections. Those radio galaxies where only a maximum 
fit of the stellar SED was possible have upper limits, indicated by arrows. The dashed lines 
represent the luminosities of elliptical galaxies with Zform = 10 taken from the PEGASE.2 
models and normalized to IQ^^Mq and 1O^^M0. Crosses mark the stellar luminosity of sub- 
millimetre galaxies from Borys et al. (2005), re-derived in the same fashion as our HzRG 
sample. 
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Fig. 10. — Rest 5 //m monochromatic luminosity versus redshift. Filled symbols indicate 

sources observed with MIPS; open symbols indicate sources without MIPS observations but 
observed with IRS. Dotted line represents the divide between mid-IR "strong" and mid-IR 
"weak" AGN, as defined by Ogle et al. (2006). 
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Fig. 11. — Ratio of rest-frame stellar near-IR (//—band) luminosity to total mid-IR (5 /im) 
luminosity against redshift (top) and 3 GHz luminosity (bottom). This ratio appears to 
strongly correlate with both redshift and radio luminosity, although both correlations are 
probably due to the loose correlation of the 5 /im luminosity with both redshift and radio 
luminosity. 
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Fig. 13. — Stellar mass in solar units plotted against redshift for the 69 Spitzer HzRGs,. 
Symbols are as in Fig. [9] 
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Fig. 14. — Stellar mass against rest-frame 3 GHz luminosity. There is evidence of a slight 
trend, such that the most massive galaxies host more powerful radio sources, though, ac- 
counting for the large number of upper limits in the plot, a formal survival analysis shows 
only a 50% chance of the two obscrvables being correlated. Deeper data and better coverage 
with the longer wavelength Spitzer cameras would show conclusively whether or not this 
correlation exists. 



